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I. INTRODUCTION 


The advent of the digital computer as standard 
equipment cn toard virtually all modern Naval ships has 
opened the field of Digital Computer Contrcel in almost all 
aspects cf ship life. The computer has been a viable asset 
in fire ccrtrol systems for years and has been used 
extensively for aids to ship maneuvering control in the form 
es NTDS (Naval Tactical Data System) readouts. The 
declining ccsts cf generai and special purpose computers has 
made their Loception as a Manpower replacemert or 
augmentaticn a reality. Their high speed and accuracy can 
Make them fperform functions with much greater safety than 
previously attainakle with time proven fand sometimes time 


weary) “seaman's eye." 


This then is the basis for this thesis; a study of the 
types of maneuvers that can be handled more accurately and 


safely than pfresentiy being accomplished. 


One area of study is the total Replenishment At Sea 
(RAS) protlem including both course and speed control for 
tke apprcach and alongside phases. This situaticn has 
always been cne of extreme danger due tc the collision 
Fotentials invclved. However, other dangers are jnvclved in 
the on deck evclutions when the ships are not kept at a 
fairly cecnstant distance. Sudden violent maneuvers may 
cause the stress on the lines between shifs to increase 
enough tc cause the lines to part. The reality of this 
danger is readily apparent to anyone who has ever seen a 
Marila line cr steel cable part or a kingpost shackle break 
or a kingfost suddenly rend under these extreme stresses. A 


System which will minimize these potential dangers is well 


17 


werth investigation. 


Of course with a digital computer, the algorithm for 
RAS can Ee modified or replaced to enable its use as a 
Maneuvering contrcl device for other situations such as 
formation steaming or single ship navigation transit 


eentrol. 
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IT. MODELING 


A. MARINER CYNAMICS 


In the ccnception of this thesis, realistic models of 
modern destrcyer hull configurations were sought. This 
search proved fruitless. The hydrodynamic coefficients for 
present cay destroyers are not currently available. 
However, scme naval and civilian research is presently being 


conducted to cktain these coefficients. 


A cecmplete set of these coefficients is necessary for 
any maneuvering control system design. A hull configuration 
which has teen under continual study with well defined and 
verified hydrodynamic coefficients was chosenf{!]j., This hull 


is commonly referred to as the "Mariner" hull. 


The development of the equations of motion in six and 
three degrees of freedom have been well documented{*]. The 
mcdel used fer this thesis is the eguations of moticn in 
three degrees of freedom linearized with second order and 
higher terms eliminated. These equaticns are characterized 
by dependency on small perturbations about a specific 
Operating point. The maneuvers experienced in the fcllowing 
charters dc not entirely meet this criterion. The 
inadeguacy and shcrtcomings of this model are of little 
ccnsegquence Lecause relevant hydrodynamic cacefficients are 
not available, and the methods presented can be applied to 
any ship tyre. 


The development of the model is readily available to 


the interested reader [3]; only a summary of the equations 


US, 


and their ccrresfonding hydrodynamic GEctricienes are 


presented kere. 
The equations of motion used are as follows: 


X.-&)U+X (U-u )+X d=0 

( u ‘. 9? d 

(Yo-u) V+Y V#(Y -m) R+Y.R+Y d=0 
V v Ee rc d 


(No-I )F4N R4#N.V+N V+N -d=0 
ie Zz 2 v Vv d 


The direction and sense of the terms in the akove 
equations are shown in figure II-1. Letting: 


a =E-Y. 
11 Vv 
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.) 
RE eS Ee Ee ee ee ee ee ee 
Se Se A i ip gl 


> 


Figure II-l 


Direction and Sense of Hydrodynamic Terms 


Ze 


The equaticns cf motion can Le written as: 


a Atb Atc Ata Btb Btc B=IF1 
Meee ii 21060:«d1tC2 


ft te 
a Atb Atc Ata Beb Btc B=IF2 
12 1z nz Boe a2 2 a2 


a Ctb Cte C=IF3 
33 33 


33 
or: 
nN tt 
a Ata_€E=11 
<a 
rT tt 
a Ata E=I2 
Z 
t 
a C=I3 
She 
where: 
IF1=-Y <¢d 
IF2= ke 
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ie3e- f Tad 
u 


and: 


I1=-k &-c A-b B-c BtIF1 
1) 1 ee ae; 


I2=-t ~B-c A-b B-c B+tIF2 
Dow 22 


By sclving this system of equations, the following 


relaticnshifs are established: 


A=(a_a_TIiw-a_.a_12)/DEL 


ey 35 3.21 


B=(a a I2-a a 11)/DEL 
Ad 33 Sac AZ 


DEL=é4 (aad. -a a ) 
Sa eit 22 12 21 


which yields the solution: 


- tC. it 
V=A=V “f A dt 
0S, 
t. t t T | 
+ 


B dt=0 fo cB (oy +f Edt jdt 
te 0 to ty 


q € ois 
U=C=u f E che 
0+ 


The space cocrdinate system is defined as fcllows: 
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Y=UeSING+Vecosyp 


X=U «COSP-VeSIND 


where: 


t 
X=X +f X dt 
¢) te 
t 
Y=Y f 
Of, 


Table II-7 summarizes the symbols and nemenclature used 


dt 


> a) 


in the fcregcing abbreviated solution of motion in three 
degrees or freedom. The applicable characteristics of the 
Mariner hull are presented in table II-2 with the 
corresponding nondimensionalized hydrodynamic coefficients 
and the DSL computer program variable names delineated in 
takle II-3. | 


Computer Program #1 is the basic DSL program that was 
developed frem these equations of motion. This prograkz uses 
two ships to illustrate the turning characteristics cf the 
Mariner hull fcr various rudder commands. Figure [I-2 shows 
the difference between a step mcdel and a ramp nodel rudder 
in a gecgraphic flot. Figure II-3 shows the correspcnding 
difference in yaw. 
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TABLE II-1 


SYMBOLS AND NOMENCLATURE 


derivative 


respect to 


derivative 


respect to 


derivative 


transverse 


derivative 


transverse 


derivative 


tc angular 


derivative 


to angular 


derivative 


Ep =e qe SS SSS 


CreLOngituda nal trorce MITER 


longitudinal acceleration JU 


of longitudinal force 


longitudinal velocity 


of lateral force 


acceleration V 


of lateral force 
velocity V 


of lateral force 


acceleration R 


of lateral force 


velocity R 


cf lateral force 


to rudder angle d 


with 


with 


with 


with 


with 


with 
U 


respect to 


respect to 


respect 


respect 


respect 


derivative of yaw moment with respect to 


transverse 


acceleration V 


derivative of yaw moment with respect 


tc transverse velocity V 


ZS 


roe 


«je 


TABLE II-1 (cont.) 


SYMBOLS AND NOMENCLATURE 


Definition 


derivative of yaw moment with respect to 


angular acceleration R 


derivative of yaw moment with respect tc 
angular velocity R 


derivative of yaw moment with respect to 
rudder angle d 


yaw angle acceleration 
yaw angle velocity 


initial velocity of origin cf body axes 


relative to fluid 


transverse acceleration of ship axes 
relative to fluid 


transverse velocity of orign of ship axes 
relative to fluid 


hydrodynamic longitudinal fcrce 


hydrodynamic lateral force 


longitudinal acceleration cf Ship axes 
relative to fluid 
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TABLES li-t (cont .} 


SYMBOLS AND NOMENCLATURE 


symbol Definition 
U 


aoa velocity of ship axes relative 


yaw angle 
actual time 
nondimensionalized time 


longitudinal distance that the ship center 


of gravity is forward of the ships axes 


longitudinal velocity of ship axes relative 


wenptluids (Cperating pcint) 


PAB LE £1=2 


CHa AGLERT SLES OF MARENER-TYPE STUDY SHIP 


Length, ft. 5. 
Beam, £t 76.0 
Draft, ft 29.75 
Displacement, tons 16,500. 
Block coefficient, S 0:6 
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TAB ioe tas 


NCNDIMENSICNAL HYDRODYNAMIC COEFFICIENTS 


Computer Program Nondimensional 
Coeficient Variable Name Value 
(X.-b) MXUD -0.0085 
u 
X XU -0.0012 
u 
iy vay -0.01243 
v 
(v2) MYVD SU SS) 
Vv 
(Y -mu ) MYR -~0.0051 
(Y.~-ex ) YRD =O. 0100127 
g 
3 YDELR +0.0027 
N NV “0.00351 
v 
N. NVD -~0.000197 
v 
(N -mx u ) NR -0.00227 
(N.-I ) IZNRD -0.00068 
ee SZ, 
Me NDELR -~0.00126 
X Sa -~0.0000462 
n 
iy a ~0.0000052 
n 
N --- +0.0000026 
n 
X XDELR 0.0 
d 
NOTE: x =C€.0 
g 
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NCNCIMENSIONAL HYDRODYNAMIC COEFFICIENTS 


Values Eased on the following operating point: 
eg? (15 Kts) 


2g 


»*O§ OC 'Ce ct "0g C005 vo OL Go°oeg Go ‘0? Cott ou 'O 


(2) x) BINvWISIQ 


A3zLINVWHXON 


sS 
x 
Sadan awyy © ysogndg AILS 
x 


uosTtazeduog 
oTyderboay Tapow duey pue days Taeppny 
c-ILl eanbty 2 


30 


ZOSTL C4 -wwQ Qionkse zyw 


00°08 02 °O¢ satay N0°O3 MOO: 30°Oc 
Le ee: | Ee a ae a Parr Sra 
(Or 93S) Bx i 
i ee 
0 eee 


SUT, “SA so.UueTeFIIG 
MeK Tepow duey pue days rzeppny 
€-II sanbtq 


cO°Or 


mM 


Og °T- 


HIF A-uuweuzyu(awy) 


a 


Bae RUDDER RESPONSE 


The previous section indicates a marked difference in 
behavicr retween step and ramp rudder models. This prempted 
an investigation into realistic rudder modeling which would 
fulfill tte reguirements of limit stops and maximum rcudder 
rate. . 

NSRCOC{*]} has modeled the rudder of the DD-931 Class 
Destroyer. The basics of this model are presented in the 
block diagram of figure II-4. 


The first limiter models the rudder stcfs which fcr the 
Mariner are +30 degrees. The second limiter models the 
proportional band of a variable-displacement pump by 
ligiting its maximum percent stroke. The limits for this 


nonlinear element have been found to be +7 degrees. 
The transfer function (K ys) acceptS an input error 
9g 


Signal of up to 7 degrees, converts it to a rudder rate, and 


integrates the rudder rate to obtain rudder angle. Letting: 
d = Maximum rudder rate (2.0 degrees/sec) 
D 


d $= Maximum error input (7.0 degrees) 
emax 


The system gain can Le defined as: 


K =a, 
g - emax 


2507 120 
Q.2€5714 /sec 


5 


werbetq yoo,Td Aeppny 


P-II oeazanbty 
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Tc convert this model to the required nondimensionalized 


form, the fcllcwing manipulation is required: 


K §'§ = K eL/u 
g g 1 
= 0.285714e527.8/(15¢1.689) 
= 5.95224 
wheres; 


L = ship length 


a, = cpéerating point speed (15 Ktse1.689 ft/sec/Kt) 


Computer Program #2 is the DSL program which grodels 
this system. The curves of figure II-5 and II-6 exhibit the 
responses cf various step rudder commands. These are 
tabulated and cress referenced in table II-4. 


These responses show the characteristics of a realistic 
rudder in that the rudder is never allowed tc slam intc the 
stcps. They exhibit the time delay between command and 
response which is a function of the rate of response (2.0 
degrees/sec). A centrol system design with this scheme is a 
much more difficult problem than one with an idealized 
rudder (step resfonse) because the entire rudder ccntrol 


system beccmes quite nonlinear. 
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TABLE II-4 


RUDDER COMMAND AND RESPONSE 


Rudder septal ise 
Figure Curve Command (deg) Condition(deqg) 


ss 1 +30.0 -30.0 
rs a2 +25.0 SA 
it ='> 3 +20.0 = P4100, 
fi 5 4 +15.0 = 7 (8, 
1S) 5 +10.0 are O 
i= Ss 6 Poi), = 5-0 
il —6 1 +3C.0 O50 
= 6 2 25.0 0.0 
tel = 6 3 +20.0 0.0 
ET=6 4 riae@ 0.0 
TI-6 5) +#10.0 0.0 
ri © 6 + 5.0 0.0 
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Figure II-7 pertrays a complex model of a gas turbine 
propulsicn flant{>j[*]j. This model cantains the elements 
required for a complete dynamic study of the systen. For 
the purpese ci this thesis, such a complicated model is not 
required if the overall input-output relationship can be 
estaklisted. 


Reference 5 establishes an . output speed (0) 
relaticnship for a step input of desired speed (Ud) and is 
redrawn as figure II-8. The relationshifp aprears to be that 
of a first or seccnd order system with a time delay. 


The system equations for a first order approxization 
With a tine delay say be written as: 


SPLIIN = K*SPDDES*¥e 

SEDERR = (SEDIN-SPEDCUT) *G 

SELOUT = f severe at 
Which yeilds the transfer function: 


i) 
SELGUT(s) G*K*e 


SPLIIN(s) $+G 


_ Which is bleck diagrammed in figure II-9. 


SPODES io e _ SPDouT 


| Figure II-9 
Propulsion Plant Low Order Model Block Diagram 


Sits) 


| (2 SP 2) + 


(“Pry ave'v Sa) Fo 7 
t 


werbetq yooTd yueTd uotstTndorg sutqany, ses 


L-II oanbty 


s/t 


yy) 
(vp f+ (4y -2L) = 21 
+ Jif 
D1 9-1): AL 


(,92% 4, dA), Gd AD = 91 
d:g-D=dO 


ay) 


O'0ce 


O'O82 


(°93S) Fwit 


asuodssy ders Tepow zepzAO yhtH AUeTA uoTsTndorg 
8-II eanbt 4 


O'or?2 


O°004¢ 


0'09/ 


O'O? ! 


0°08 


O’Or 00 


bo? 


SoZ. 


Cie 


Viz 


NALOODSE (x kw) 


40 


From figure II-8, the time delay, system gain, and time 
ccnstant can be estimated as: 
T = 4.88 sec 
0.9877 
0.092 


QQ 
ll 


Computer Program #3 was used to okttain the ster input 
response. The original complex system output and the low 
order approximaticn are compared in figure II-10. As 
indicated in this figure, the two responses are very close. 
Considering the linearized approximations made in the 
eguations cf mction, this response 1s accurate encugh for 
system study use and is used as the model fcr speed ccntrol 
in chapter III. 


Similar methods may be used to oftain simplified 
low-order acdels for other. high-order propulsion systems now 
in use (e.g. pressure fired boiler systems, 1200 1b. 
systems, etc.). They may not, however, simplify to a first 
order afrreximaticn suitable for system study. A method of 
computer determination of low-order models of high-crder 
Systems is ccntained in ref. 7 and may be mechanized for 
this purpese. 
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D. EXTERNAL FORCES 


The modeling of ship dynamics cannct be complete 
without the introduction of external forces which perturb 
its respenses. These forces are caused by many factors and 
some are mcre relevant than others in the scope of this 
thesis. The two that are considered can cause substantial 
perturbations that must be modeled and eventually acccunted 


for in tke ccatrcl system design. 


Whenever two ships operate in close proximity (less 
.than 250 feet), suction and pressure forces between hulls 
are present. Studies have peen conducted on the Mariner 
hulif[!] which have produced data for ccnstruction of a 
fanily of curves fcr two ships passing on the same heading. 
No data has teen found for the cases of two ships not cn the 
same heading. Other restrictions on the work presented in 
ref. 1 are that the ships are of the same tyfe and of 
Similar hull ratios. 


Interactive effects vary as the square of speed. 
However, this is only true if both ships are at the same 
speeds. The interaction modification factor is based cn the 
normalizec speed of 15 kts. This factor can thus be written 
acs 


2 
SPDE = CDOT! 


Exact effects on the interacticn forces and moments in the 
Situation where the ship's speeds differ are not available. 
This is inconvenient since during the approach phase, the 


nciImalized speed of the approach ship (ship 5B) can fe as 
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high as 1.5. If the effect on ship B is as stated above but 
With its cwn speed causing the interaction modification, the 
interactive forces and moments can be 2.25 times greater 


than withcut speed considerations. 


Withcut the ability to pin down this relationship, it 
was ignored in the development cf the contrcl laws presented 
in this thesis. Appendix C was written with the expressed 
intent of illustrating the effect ope modifying the 
interactive forces and moments to the extremum case 
menticned akcve. It must be realized that this case is not 
ccnsidered likely in that it is felt that the interactive 
forces and mcméents modification cn ship B are more apt to be 
caused Ey the speed of ship A. If this 1S so, since ship A 
is kept at a constant 15 kts., the interactions need ano 
mcdificaticn fcr speed consideration in this thesis. 


Reference 1 also gives a method of modifyince the 
interactive forces and moments based on different ship 
lengths. FCr ease of computation and graphical 
presentation, the two ships were considered of equal 
lengths. Tc modify this to ships of dissimilar lengths, the 
resulting hydrcedynamic derivatives must be mcdified as shown 


in ref. 1 (alse shown in appendix Cc). 


Since neo closed form expression existed for these 
forces, tke fagily of curves reproduced as figures II-11 and 
II-12 were guantized in subroutine SLOPES (an adaptation of 
the subroutine of the same name from ref. lieu (A DPend1x .9c 
contains a curve fitted subroutine that was compiled after 
tke completicn of the research on this thesis. It was not 
used for any design qc simulation runs except for those 
presented in that same appendix.) An interpolation 
algorithm is used to approximate the intermediate values 
between cuantized values and between the curves of the 


family. 


au 


Lateral Distance DY (Ft) 


-400 -200 0 200 400 


Note: To convert to normalized lateral distance DY - 
Divide by the ship length L. 


Figure II-ll 


Family of Interactive Y force Curves 
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* 
Lateral Distance DY(Ft} 


-400 -200 Q 200, 400 


* 
Note: To convert to normalized lateral distance DY - 
Divide by the ship length L. 


Silefcbelsy JE 


Family of Interactive N Moment Curves 
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The main purpcse of this subroutine is to compute the 
interactive forces between ships in the replenishment at sea 
Situation anc output the values for perturbation of the 
ccentrol ship only. Even though both ships are affected by 
‘these perturkaticns, a one ship control system which is 
effective regardless of the other ships moticn is 
considered. Ccnseguently, the interactive forces on the 


seccnd ship are ignored. 


Subroutine SLCPES is contained in appendix A. Figure 
IiI-13 is a geographic plot of the two ships rfassing at 105.6 
feet with their rudders amidships (0 degrees). Figure II-14 
and figure II~15 show the magnitude of the lateral force (Y 
ferce) and rctational moment (N moment) of the reference 
Ship on the ship making its approach (control ship). The 
reference ship is at 15 kts. and the apprceach ship is at 
Pees KtS. The control ship starts its approach 5 ship 
dengths (2639.0 feet) astern and 0.4 ship lengths (211.12 
feet) lateraily displaced. The most graphic portrait cf the 
effects cr these fcrces and moments appears in the yaw 
Changes which are presented in figure II-16. From these 
figures it beccmes readily apparent that these perturbations 
cause viclert motions of the ship which must be acccunted 
for in any ceontrcl system development. Throughout the 
development of stch a control system in chapter III, these 
forces and mscments are considered inherent in the model for 
Hee GCOtCrcl. 
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2. Waves 


The nodeling of sea state in the form cf waves and wave 
interacticns has cccupied the time of many researchers 
cot ? j. The exact formulaticn of waves will act be 
accomplished in this thesis. Since the main concern here 
Will be toc test the control scheme developed in chapter III, 
a much sigzplified wave generator can be used. To introduce 
the required experimental perturbations cn the designed 
control system a periodic wave system with a fundamental 
frequency ard its second harmonic is used. Some small 
Candom wave properties are introduced that ride on these two 
sinusoids. A simple expression of this combined wave can be 


written as: 


2 
W = WFesin(WE)+(ELeWF /WL) esin(2¢WE) +WFeWRVesin (WE) 


iW denctes the Wave 

WF ais the Wave Force 

WE is the Wave Encounter radian Frequency 
WL is the Wave Length 

WRV is the Wave Random Variable 

PI is 2=.1415926 


With this wave as a basis, a method cf acdeling this in 
the dynamic environment of the total ship simulation was 
defined. The modeling includes the intrcduction of this 
wave intc the three degree of freedom eguaticns of mction. 
To accomplish this a set of defining relationships were 
develcped. First the general wave directicn is input to 
establish the direction of the wave encounter on the ship. 
If the ship direction is YAWDP2 and the wave directicn is 


HD, the expected wave direction is defined as: 


EWC = WD-YAWDP2 
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Next the wave enccunter frequency (radian frequency) can be. 
established with knowledge of the ships ncrmalized true 
speed (CDCT2), wave length (WL) and normalized wave velocity 


(WV). The wave encounter frequency (WEF) is then: 
WEF = zePIe (CDOT2+WVecos (EWD))/WL 


The total wave encounter (WE) is ncthing more than the 
wave enccunter frequency (WEF) times time. This gives the 
wave enccunter radian frequency required in the simple 


expressicn fcr the combined wave previcusly shown. 


This dces not ccmplete the task, since the individual 
wave forces of each degree of freedom must te derived for 
this general wave expression, namely the ccmponents of WF. 
Again a much Simplified version of the more complex real 
life wave forces were used. The X and Y forces are 
considered first. These can ke modeled as cosine and Sine 


functicns cf the expected wave direction (EWC) such that: 


WE X 
WEY 


WFecos (EWD) 
WPesin (EWD) 


where WF is the total wave force of the enccuntered wave. 


The rotational N forces are a little more difficult to 
establish. Ey ccnsidering that no rotational forces are 
created Ey a wave directly on the bow or stern, or directly 
cff the teagz, and that it is maximum when the wave is at 45 
degrees off the bow Or stern, a much Simplified 
apEroximaticn is developed. Realizing that this method is 


very crude, the N force can be written as: 


WEN = WFesin (2eEWD) 


To add mere creditability to the wave mcdel, a random 


ug 


variable is added to the wave force at the waves fundamental 
freguency. A gaussian (normal) distributicn was chosen with 
a zero mean and a standard deviation of one-tenth the 
Maximum allowarle force of WF. A zero mean signifies that 
the expected amplitude of the random wave is 0.0, while the 
standard deviaticn signifies that 68% of the random waves 
will have agplitudes less than one-tenth of the maximun 
allcwable ferce of WF. Also, 94% will have amplitudes less 
than one-half the maximum allowakle force of WF. This small 
added perturkation allows for verification of the model 
Simulaticn wath a stochastic force, which in turn adds 


creditarility to the developed ccntrol systems. 


What remains is to define the total wave force (WF). 
It is impcertant not to fall into a common Simulation pitfall 
Which inevitably causes unneeded design changes. A sea 
state does nct increase at an infinite rate. It therefore 
is incorrect to start a simulation with initial conditions 
set for calm séa and immediately introduce a high sea state 
perturbaticn. The initial large perturbaticn transient can 
give results that are not only unrealistic, but can cause 
the model and control system to produ3e unstable results. 
This is especially true in this case since the linearized 
(stall perturkaticns about an operating pceint) equaticns of 


mcticn are uwsed. 


With this in mind, a ramp feed in of the wave force 
With a limiter at the desired maximum wave ferce (WFMA) was 
used. The slcpe of the ramp was established to impart 
Rinimum initial transients, yet increase the wave force to 
an accertable testing level within the time frames cf the 
Simulaticns cf chapter III. The slope is designed such that 
the maximtum wave force is reached in 94.815 seconds actual 


Simulaticn time (4.548 seconds problem time). 


Computer program #4 was used to verify the wave action 


50 


model. Table II-5 on page 56 indicates the figures freduced 
and changes in input wave length (WL) and wave direction 
(WD) for each run. The input parameters that were constant 
for all runs are tabulated below: 

YARKDE2 = 0.0 

CICT2 = 1.5 


* 
WS = 5.0 
WFMA = 0.1137 


x 
NETE: WS is the unnormalized wave speed. Conversicn to 
ncrgalized wave velccity is: 
WY = WSY15.0 


Intrcduction of the wave forces is accomplished by 
multiplying the established wave forces by the rudder 
hydrodynasic coefficients for the individual reference 
directions. This effectively scales the wave forces tc the 
ship model tckéing used. The wave rfrorce result 1s coded in 


tke ship simulation program as follows: 


IF12 = KA1eD2+YY2+KAIOWY 
IF2Z = KELODZ+YN2+KBIOWN 
IF32 = KC1eD2+NC2+KC1 oWX 


Detailed results of the wave farce effects are given in 
Chapter III and will not be dealt with here. 
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* 
NOTE: WL is given in ship lengths 
WD iS given in degrees 
Curve numbers of all runs corresponding to 
wave force components are: 


Force Curve 
WX Be 
WY 2 
WN S 
Table II-5 


Wave Simulation Listing 
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ite CEPLENTSHMENT AT SEA 


Aj. HEADING CCNTROL 


tes Control Chcice 


Many studies involving replenishment at sea (RAS) have 
treated the preblem asa multivariable system{ 19 j[11)[12] 
{43}. Academically, there is nething wreng with this 
approach. However, aS a practical system it leaves much to 
be desired. The key drawback in the multivariable system is 
the inescapakle dependency cn a command and control link 
between tke replenishment ships. The unreliability of UHF 
communications at these close distances is a much 
experienced phenonencom. It is felt that any knowledgeable 
ccamanding cfficer would not entrust the safety of his ship 
to such a guestionable link. An alternate method which is 
described kLere is a modern extension of the long trusted 
"seaman's eye" ccncept, where the senscrs and ccntrol 
devices mtst be self contained on the Ship making the 
@meecach (tereafter referred tc as the receiving ship or 
Sap 8B). 


In all fresent day RAS operations, the ship on which 
the apprcach is being made (hereafter referred to as the 
supplying shir or ship A), must maintain the replerishment 
course and speed. The receiving ship accomplishes the 
Maneuvers tc Baintain station relative tc the supflying 


She 


The parameters which are presently measured "visually" 
are relative position (in both the X and the Y directions), 
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relative heac (usually in reference to crdered replenishment 
course), and relative motion in the X direction (fcr speed 
matching). ‘These parameters are usually visualized by the 
ccnning cfficer who in turn gives corrective orders to the 
helzsman. The helusman must then translate these verbal 
orders into rudder and speed commands through the helm and 
lee helm ccnscles. The accuracy of the execution of the 
conning cfficer‘s ocrders 1S extremely dependent cn the 
ability cf the helmsman and throttleman. This system can be 
guite effective, and it can also be quite disastrous. MThis 
fast reacting and constantly changing environment lends 
itself tc Ereakdcwn in communications and manifests the 
inability of scme individuals to cope with the required 


critical ran-sachine interfaces. 


To eliminate these problems, present state of the art 
digital cecmruters and sensors are available for immediate 
itplementaticn of a completely automatic ship Gentnor 
system. Such a control system may be installed on 
individual ships and be used for RAS without the reguirement 
of having the matching installation on the cther ship cf the 


heckup (another drawback of the multivariable approach). 


One cf the many pitfalls that may be encountered in 
digital simulation is the reality of the parameters that are 
measurable ir the real world situation as opposed to those 
that are incidently available in the simulation. With this 
fact as a keynote, Subroutine RBMEAS (Range and #é¢aring 
MEASurement) was developed. This subroutine, as listed in 
appendix A, defines the forward (FWD) and after (AFT) 
relative and true bearings, and ranges from the receiving 
Ship to similar points on the supplying ship. Figure III-1 


delineates tke terms used in the Subroutine. The SDFn terms 
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Measurement Techniques 
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are used tc position the bow and stern sensors and 
reflectors cn ship B and A, respectively, in geographic 
cocrdinates as a function of ships head. The FWD distance 
on the xX cccrdinate is ADPX and the Y cocrdinate is ADFY. 
Sigilarly, tte AFT distances are ADAX and ADAY. R1 and R2 
are the FWE and AFT ranges measured by a highly accurate 
rLanging device installed on ship 8B. This same tLCanging 
device, if properly provided with a pinpcint reflector on 
the supply ~ship(ship A), will give accurate relative 
bearings FWD and AFT., B1 and B2 respectively. The distance 
between senscrs may be varied, but as a rule should te kept 
as far apart as pessible to allow maximum sensitivity. The 
distance used in this thesis is 1.0 (one ship length), and 
the distances were considered the same for beth ships. This 
is net a necéssary condition and may be changed to suit the 


Ssituaticn. 


Subreutine RBKEAS assumes highly accurate sensors in 
both rance and Fearing measuring ability. Such sensors are 
presently available in the form cf Radar altimeters{1*] and 
Laser ranging devices. Another possibility for a measuring 
method is a single sensor time sharing to cbhtain range and 
bearing tc kcth reflectors from a single device. Such a 


Sirgle senscr scheme is sketched in figure III-2. 


Once tke FWD and AFT parameters are available, they may 
then be used to determine other desired quantities. 
Subroutine HDGRAS ({HeaDinG control for RAS) was develored to 
output tke desired heading corrected for heading difference 
Gfeeeship A and B and the prejected correcticn for distance 
errcr. This subroutine is listed in appendix A. The center 
Cange and tearing are the average of the FHD and AFT range 
and bearing cutput from Subroutine RBMEAS. 


The additional heading due to distance is projected as 


if ship E saintained its present course until it was 
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perpendicular to the center of ship A. The reasoning Eehind 
this is illustrated in figure III-3. If the present course 
will cause ship B to arrive on the station desired (DS), no 
heading chance 1S reguired. The expression for PSIACC (J 
Additioral heading due to Distance Correction) is: 


| pstapc = RSENS (DDC+DAeSIN (AA1) ) 


where: 
RSENS = Range SENSitivity gain 
DDC = Distance Desired Corrected for side of apprcach 
DA = center Distance Absolute (range) 


AA1 = 3€0 degrees - relative bearing of center position 


The beadince difference of ship A and B is desired 
Since, even if the range when alongside is correct, a large 
disparity in héading cannot be tolerated. It is realized 
that some heading difference (crabbing) is necessary to 
Maintain tke distance. This crabbing is due entirely to the 
pressure ferces modeled in chapter II. This heading 
difference is found by computing the difference in the 
perpendicular prcejection between the FWD and aE 
measurements and finding the arcsin of this difference 
divided Ey the distance between sensors. Figure IIiI-4 


indicates a sample of this procedure. 


The expressicn for total desired heading is given as 
Eclicws: 


PSIALC+NTSENSePSIDIFtPSIE 


PSIDES 
where: 


PSIDIF = Y additional heading due to heading DIFference 


WISENS WeighTed heading difference SENSitivity qain 
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Figure III-4 


Heading Difference Calculation 
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PSIDES = § (heading) DESired 


Thretghcut the subroutines and main DSL programs, the 
Function CEGRAD (conversion of degrees to radians and 
radians tc degrees) is used freely. An explanation and 


listing cf this functicn are presented in aprendix A. 


“The angular velocity of the receiving ship's head is 
alsc cf cencern in the RAS Situation. This guantity may be 
thought of as Similar to tachometer feedback in a Simple 
serve ccntrcel system; and is necessary to damp out the 
response (the responses of this control system without this 
feedback igs presented in the latter section of this 


chapter). 


The desired rudder command is a combination cf the 
desired beading, angular velocity feedback, and a rudder 
gain as fcllcws: 


Desired Rudder = (YAWD2-PSIDEL+BDOTFB) «GN 
where: 

YAWDZ = heading of ship #8 (in degrees) 

PSILED =ESIDES (in degrees) 

BDOTFB = VFBGeBLOTZD 

VFG = Velocity FeedBack Gain 


BLOTZD = angular velocity of ship B heading angle (in 


degrees/sec) 
EBDOTFB = angular velocity FeedBack 


RGN = Rudder Gain 


The convention for rudder response dictates negative 


Zl 


rudder as kéing right rudder, which causes positive yaw. 


This necessitates making the desired rudder the negative of 
the forcing function and feedback quantities. The block 
diagram cf figure III-5 presents the control loop 
measurement inputs to desired rudder command. 


from 
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Optimization 


Thus far the control chcice has identified four gains 
(RSENS, &ISENS, RGN, VFBG) that must be set for [freper 
position attainment. The necnlinear nature of the systen 
which afrears in the form of distance measurements, 
interactive forces and rudder modeling do not alicw for 
straight forward determination’ of these gains with normal 
optimal ccntrcl theory. , 


Grossly nonlinear systems require special handling to 
determine prefer gain settings. The method chosen for this 
purpose was an optimization algorithm develored by M. J. Box 
(ercgrammed locally as subroutine POXPLX). fTfhis subroutine, 
listed in appendix A, was used to locate the cost function 
saddle pcint in four dimensional space (the dimensicns being 
the previcusly mentioned gains). The drawback asscciated 
with this wethod is the necessity of iterating the complete 
ncnilinear sinulaticn within function FE for every evaluation 
of the ccst function. The gains sought were found, but 
unfortunately only after 2 1/2 to 3 nours of CPU time with 


every 400 iterations allowed. 


The weéechanics involved in optimizing the chosen cost 
function include required sub-calculations in many functions 
and subkroutines. Figure III-6 is a flow chart which 


demonstrates the steps, subroutines and functions required. 


Initial optimization was accomplished for one set of 
Mier ai conditions. By looking at the RAS situaticn, a 
rrckable set cf circumstances were envisioned. The scenario 


setting is the approach phase where the rerflenishing ship 
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Optimization Flow Chart 


TS 


starts a wide approach at 0.4 ship lengths (211 feet) 
lateral displacement and 5 ship lengths (2639.0 feet) astern 
of the suprly ship. The desired final position is alongside 
at 0.2.ship lengths (105.6 feet) lateral separation. Eoth 
Ships have tke same initial heading (YAW angle). The supply 
ship is at 15 kts. (1.0 normalized speed) and the receiving 
Ship makes its approach at 22.5 kts (speed control will be 
covered later in this chapter). 


Normal ccstinc¢c of displacement error with the integral 
of the squared error (ISE) was considered as the 
optimization tcol in subroutine BOXPLX. However, this type 
cf perfctmance measure would weigh the initial displacement 
errer equally with the final position error. This problem 
can be circumvented by comparing the displacement error to a 
pre-ccmputed reference track instead of to the desired 
displacement. For the envisioned scenario, it was conceived 
that the ccst function should weight the distance 
displacement teavier when the ships are alongside than when 


the apprcach is started 5 ship lengths astern. 


This wés accomplished by using the integral cf time 
times the aktsclute error (ITAE) as the optimization 
performance measure. The reference displacement was 
considered tke desired position displacement. rue Object 


furcticn can then Le written as: 


te 
CBd -/ te|DD<-ADY] dt 


t, 
where: 
DD = Desired Distance 
ACY = Actual Displacement in the Y directicn 


t = time 
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A perfcrmance measure that is designed to obtain good 
performance must also take into account other factors 
besides just fosition accuracy. Consequently, another cost 
criterion was decided upon which would also set the gains to 
reduce tte agrcunt cf rudder contrcl required when alongside. 
This particular feature is derived from the desire net to 
over contrcl with the rudder in such clcse proximity to 
ancther vessel. The inclusion of this term in the 
performance measure is weighted by unity while the distance 
accuracy 1S weighted by a factor of 10.0. This will tend to 
allow rudder action if the desired pesition is not 
Maintained. The final approach phase cost function for 


oktaining cptimum gains has the form: 


t¢ 
CEJ = i te(10.Q0e«{DD-ADY|{[+1.Q0¢e{D2]) dt 
te. 
where the adéiticnal term is: 


D2 = rudcer response of the replenishing shir 


Cc. Resuits 


In the process of deciding on the best gain definitions 
previously mentioned, many optimization runs were made. 
Each set cf gains were then simulated in a corresponding DSL 
program to cktain performance confirmaticn. Many of these 
runs did not live up to expectations; causing re-evalvation 
of the ccntrcl scheme until the one presented in this thesis 


was formulated. 


Table III-1 shows the input upper and lower limits of 
s€arch (EU, EL), starting value guess (XS), optimum gain 
settings (Output) and associated object function value (CBJ) 
for 20.0 seccnd normalized time simulation run in function 
Es These values were then introduced into the DSL 
Simulaticn pregram listed as program #5. The results of 
THiS Si@ulaticn are shown in figures III-7 thru III-12. The 
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gecgraphic flict cf figure III-10 indicates excellent 
positioning in the lateral direction while the rudder 
response cf figure IITI-12 shows that it settles out to a 
fairly ccnstant steady state value as the ship settles into 
its desired fcsiticn. The time coordinates in all flots are 


shown in actual full scale time. 


d. Contrcel Testing 


2S 


Now that the "proper" gain settings were cbtained, more 
extensive testing of the control system was required. Three 
different tests were contemplated: (1) allow a large 
perturbaticn turn of the reference ship (supply ship), (2) 
start approach of the receiving ship (ccntrol ship) from 
different initial conditions cf lateral and hcrizontal 
displacezents, and (3) induce external perturbations in the 


form of wave fcrces. 


The first test was Simulated by turning the reference 
Ship by ncermal rudder action of figure IifI-13. Tiaas Wee oie 
With 5 degrees rudder accounted for a total reference yaw 
Change of 15 degrees. The rudder action of the ccntrolled 
ship shcwn in figure III-14 was as expected. However, the 
distance taintainment portrayed in figure III-15 was totally 
unacceptable. The maximum excursion frem the desired 
distance cf 105.56 feet (0.2 normalizes distance) was £5.419 
feet fe. 105 normalized distance). Variances of this 


Magnitude cannot be tolerated in the RAS environment. 


Faced with this situation, the tact chosen was to 
re-evaluate tke gains for the new scenario which is called 
the turn phase. In this phase the initial conditions assume 
steady state pesitioning alongside such that the lateral 
fesition displacement (DY) is equal to the desired distance 
(105.56 feet (0.2 normalized) ] and that the hcrizontal 


Beeition displacement (DX) is 0.0 (alongside). 
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Scre initial ferturbation is introduced by assuminc the 
relative yaw angle when alongside is negligible. 


The same procedure was followed in obtaining gains that 
would optimize a chosen cost function. Figure III-6 still 
applies excerpt that function FEA is replaced by functicn FEB 
(listed in aprendix A) to simulate the new conditions. 


Cost function criteria change in this instance since 
the ships start at the desired pesition and optimally stay 
at the same relative positions. Also, the rudder response 
to such a large turning perturbation must be free to cause 
achievement Ci the desired position. Due to these 
considerations, the integral of the absolute error (IAB) 
performance measure was chosen for the optimization 


@eatrerion ane can ke written as: 


ty 
OEJ = / {ADY|{ dt 
t 


o 


Table III-~zZ shows the results of the turn phase 
optimizaticn and the comparison with the approach phase 
gains. Acain CSL simulation was performed using the turn 
phase scénario. Eiguees) E£iif-16 Sthru E£lI-21 portray the 
graphical results. The rudder response cf figure IItI-21 
indicates very sensitive response to the interactive forces 
shewn in figqtres III-17 and IfI-18. The lateral distance 
separaticn cf figure III-20 indicates excellent position 
Maintainment with maximum excursion error of only 2 feet 
(0.0038 norzalized). This minimal variaticn is well within 


that which can be tolerated in the RAS envircnment. 
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To alleviate suspicions that the response from the 
gains oktained in the approach phase could be improved by 
the gains ckEtained in the turn phase, a simulation of the 
approach fhase was accomplished with the new gains. Figure 
III-22 is tke graphic display of the effect of these gains 
on the approach phase lateral distance positioning. 


Careful analysis of the results thus far clearly 
indicate the need for an adaptive control scheme-to allow 
gain adaptaticn to meet the design specifications. A _ full 
acaptive ccntrol scheme for systems of this type is outside 
the scope cf this thesis. References 15 thru 223 are 
‘indicatiors cf some of the literature available for ptrsuit 


of a ccmpletely adaptive control systen. 


What was done here is develcpment of a Simple algcrithao 
to sense when the cenditions were adequate to switch from 
cneé set cf gains tc another. This may be done with the two 
sets of gains developed thus far. However, for the sake of 
Sigulaticn efficiency, a third set of gains was introduced. 
This third set amounts to a change of one approach gain 
(RSENS) which has previously been defined as the range 
sensitivity cain. The simulation efficiency is increased by 
decreasing tke time required fcr the approach phase to reach 
Steady state. A conseguence of this procedure is a 
Teenforcement of the need for a completely adaptive ccntrol 
scheme. 


Repeated simulation revealed that ccmmencing the turn 
(in effect switchihg gains), before a reasonable steady 
State was reached caused results similar te those shcwn in 
figure IIIJ-zzZ. An increase in RSENS to a value of 4.0 when 


the lateral separation error is less than 0.05 (normalized) 
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and greater than 0.005 (normalized) forces acceptable steady 
state in affrroximately 1/2 the time previously reguired 


usirg a single set of approach gains. 


Subroutine S&TCH (listed in appendix A) inccrpcrates 
this Simple acaptive gain schedule with a ccunter tmechanisn 
to sense wken steady state 1s reached. Further study 
indicated a need to damp the yaw oscillaticns to a greater 
extent if the yaw velocity (BLOT2D) exceeded 2.0 degrees/sec 
when the cairs are initially switched to the more sensitive 
enes of the turn phase. This is an artifical adaptive gain 
fer VFBG caused by computer time restricticns prevalent in a 
full scale computer simulation where both the approach and 
turn phases are desired. If the gain switching point is 
moved up in time, as would normally be the case in a real 


life situaticn, this damping increase would not be required. 


The results of the full scale Simulation using cogrputer 
program #6 are shown in figures III-23 thru III-34. The 
aprroach fhase fplots. of figures III-23 thru ITI-2@ show 
definite improvement over that previously shown in figures 
Zo i-/ ereu Iii—i2. Figure IIfI-27 indicates that the 
overshcot is reduced to 10.6 feet (0.02 normalized) as 
opposed to 113.9 feet (0.034 ncrmalized) that was prevalent 
in figure III-11. 


The ttrn phase plots are shown in figures III-29 thru 
fiI-34 and shcw responses very similar to those shcewn 
previously in figures III-16 thru III-21. The only 
Significart differences occur in the initial responses which 
are due to the incorrect initialization when the turn phase 
was Simulated individually. 
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men Varying initial Conditions 


The results oktained in the previous section are most 
gratifying rut actually incomplete. This system must work 
fer cther initial. conditions quite different frem those 
envisioned in the optimization scenario. The initial 
approach can realistically commence at points other tkan 5 
ship lengths astern and displaced by 0.4 ship lengths. 


By sigulating this systen with varying initial 
Feositions, the relative efficiency and worth of the ccntrol 
system can te cbserved. This was done in successive test 
runs whose iritial conditions and corresponding plot figures 
are tabkulatec in table III-3. For the sake cf brevity only 
these figures required to illustrate the relative efficiency 
of the ccntrcl system are included. The correspcnding 
initial cfptinizaticn simulation figures are listed for cross 
reference. The turn phase plots for all runs except 4 and 6 
exactly match that of the initial simulation and are not 


repeated kere. 


Runs 3, 5 énd 6 were accomplished to show that no 
ambiguities exist in the control scheme to prohibit adequate 
real life initial conditions. Run 3 simulates the situation 
mcst often encountered by this author in the RAS 
environment. This scenario starts the ccecntrol shir dead 
astern at 5.0 ship lengths and brings it alecngside at 0.2 
Shig lengths lateral separation. 


Run 5 1S a situation where the apprcaching ship is 
purposely placed out of position on the wrong side for 
approach. The control scheme adequately corrects the 
placement error and will do so for all cases of this type, 
Prcevided that there is adequate maneuvering room astern of 


the reference vessel (in this case 2.6 ship lengths was 
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experienced which gives 1.6 ship lengths bow to_ stern 


clearance). 


The purrpese cf run 6 is to frovide simulation fcr an 
approach from the opposite side again disproving any ccncern 
for ambiguity in the trignometric measurement scheme 
utilized. In all runs it must be emphasized that DD is the 
positive absclute distance desired and that IS provides the 
code flag for the desired side of approach. The system will 
work with DD set to some negative guantity; but the side of 
approach *ili reverse itself and the positicn placement will 
be correct, Eut cn the side not desired. 


Run 4 takes the desired distance in to 0.15 ship 
lengths (&€0.0 feet). This distance is usually the minimun 
desired by a prudent seaman. Again, even with this minimnun 
distance, the control system performs uf to desired 
Standards. Tke importance of this. run cannct be overlcoked. 
Performance cf the system at this extremum indicates that 
the gains utilized are correct for all expected conditions 
encountered in calm _ seas. Figures [III-59 thru III-64 
pertray tke remaining plots obtained in run 4. 


Icrmance in Sea 


h. Pertcrmance in tate 


The calm sea performance of the heading control systen 
is only fart of the system testing required. Of even 
greater ccncern is the adequacy of the control when sea 
state 1s intreduced. Section D.2. of chapter II models the 
three ccuponents of waves with two sinuscids and a small 
random impulse wave. These forces were intrcduced intc the 
tctal RAS simulation as shown in computer program #7. [In 
this program the wave length (WL) is set to cne ship iength 
and the wave direction (WD) is -015 degrees true. This 
scenario allcws for a port turn into the prevailing sé€a as 


1s commen practice in experienced RAS evolutions. BY 
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Binimizing the perturbation forces on yaw and lateral 
direction, a smocther RAS can be accomplished thus aiding 
safety and ccmfort during the actual transfer. The wave 
force maximum is taken as 0.05685. Runs were simulated 
“which used maximum wave forces in the range 0.1137 to 
0.65685, wave lengths from 0.5 to 1.5 ship lengths and wave 
directions 015 to -015 degrees off the initial replenishment 
course. The ccntrol system handled all of the perturkations 
well except fcr the cases of a wave length of 1.5. This 
length of wave with a. force of 0.05685 exceeded the ccntrol 
systems capability in that the steady state conditions were 
not met before a turn was commenced. Figure III-65 shows 
this instability in the lateral distance DY of the turn 
phase. It is felt that the modeling inadequacies of the sea 
state develcgrment of chapter II coupled with a simple 
adaptive gain scheme are the source of the problem. This 
Same phercrenon is covered in greater detail in the 
longitudinal fosition offset testing portion of the velocity 
control section of this chapter. 


Proclems of this type also manifest themselves in some 
cases when the wave force maximum (WFMA) was close tc the 
0.1137 value. If the sea state beccmes excessive, which 
this value represents, a different gain schedule or, at 


best, a mere complex adaptive gain scheme is called for. 


The flets produced by computer program #7 are presented 
aS a representative indication of the effectiveness cf the 
ccntrol system in the presence of a sea _ state. Figure 
III+~66 gives the yaw results of the appreach phase which 
indicates the effect of the wave action. The correspcending 
rudder acticn cf figure ILI-67 compensates to give the 
smooth lateral distance shown in figure III-68. The wave 
prcfile is shown in figure III-69 with curve 2 being WY and 
curve 3 -éinc WN. Curve 1 is the WX profile which was not 


used in this run but will be utilized in the speed ccntrol 
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section later in this chapter. Similar curves are portrayed 
for the turn. fhase. Figure III-70 is the yaw difference 
between the two ships (remembering that the reference ship 
is nct being perturbed by the interaction forces or the wave 
forces). Figure III-71 is the lateral distance DY 
maintained Ey the rudder response of figure III-72. The 
maximum lateral separation in the turn phase is 0.0037 ship 
lengths (1.S5 feet). The wave profile is shown in figure 
IIlI-73 with the same wave force curve sequence as the 


approach rfhase. 


As can ke seen from these flots, the ccntrol systen 
Operates very effectively in the presence of a sea state. 
Again, the develcfment of a much more complex adaptive gain 
scheme is required to allow exceptionally high sea state. 
It is felt that the control system presented in this thesis 
is adequate for most Situaticns that are encountered in the 
RAS envircnment. Only the extreme perturbations that chance 
would allow must be accounted for in a more complex adaptive 


gain scheme. 
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B. VELOCITY CONTROL 


One advantage derived from using the linearized 
equations cf motion is the decoupling of the velccity 
compenents frem the remaining equations cf motion. This. 
allows separation of the design procedures for lateral 
separaticn ccntrol and velocity ccntrol. Section A of this 
charter designed the Lateral separation ccntrol using the 
Simple speed contrcl algorithm shown in figure IITI-74. This 
control cutput was used directly as the ship's speed (CDCT2) 
in the mcdel simulation where no attempt was made to use the 
engine respense developed in chapter IE. Function SPDCTE of 
appendix A shcws the control used. 
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Figure III-74 
Non-optimum Speed Law 


Because of this decoupling assumpticn, any valid 


approach speed ccntrol can be used, if used consistantly, 
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for such a design. However, in the RAS environment, 
cctplete disasscciation is not possible. Recombination 
occurs in the interactive forces and moments which depend 
upcn the longitudinal distance as well as the lateral 
distance. Consequently, speed, which is directly 
responsikle fcr the longitudinal distance, has a direct 
relation to the lateral distance attainment and 
maintaingent. 

The remaining farts of this chapter deal with the 
development cf a viakle speed ccntrol algcrithm and the 
testing cf the designed systen. 


i. Type c£ Centrel 


Whenever two ships maneuver for replenishment at sea 
(RAS), the prime ccnsiderations are the time required for 
approach and the accuracy of position keeping lus 


conservaticn of fuel. 


The ncnlinear control law of figure III-~75 is designed 
to maintain a preselected approach speed for nininun 
approach time. The proper location of the switching point 
increases the complexity of the soluticn since the tise of 
Switching frcem this speed is determined by the dynamics of 
the nenlirnear fosition attainment loop, Once this pecsition 
1s reached, the speed controller is switched dawn toa 
linear perticn of the control law to allow control for 
perturbations about the operating position. However, small 
perturbaticns about this operating point can be tolerated 
and, in fact, are desired to allew for conservation of fuel. 
Selection of this dead zone is wholly dependent on the 
accuracy required for final fposition. Figure JII-75 
indicates a dead zone extending to +0.001 nermalized 
distance which in this case translates to +£0.53 feet. 


Systems fcr which fuel considerations are net a motivating 


to7 


factor say te designed without this part of the contrcl law 
ta allow finer tracking in the position loop. 
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Figure III-75 
Speed Control Law 
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The speed control law as explained abcve is shcwn in 
figure IITI-75 for an initial approach speed cf SPDO2 and a 
final estimated reference speed of SPDO1, with ADX bkeing the 
dynamic fosition feedback defined as _ the longitudinal 
distance between centers of the ships referenced to the 
controlled ship's heading. Analytically, the linear pertion 
of the contrcel law is written as> 


SEDCTR = -ADXe(SPDC2-SEDO1) +SPDOT 


Symetric continuation of the ccentrol law accounts for 
Operation cn both sides of the operating point. 


eke: 
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2. Cptimization 


Using this much simplified model of chapter IIT and the 
basic ccatrel law of figure III-75, the desired switching 
curve can be established. An optimization subroutine such 
aS Subrcecutine BOXPLX can be used to iteratively oktain the 
optimum switching position (Sd). for representative initial 
approach speeds., Figure III-76 is a flow chart cf the 
subroutines and functions required EOS speed ccatrol 
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Figure III-76 
Optimization Flow Chart 


Optimization. The major merit of this nenlinear contrel law 


Stems £rez the predetermination of the switching pcint for 
all possible conditions of initial speeds. This apriori 


or ee rr oe 


knowlege allows for offline computation cf the switching 
pesiticn prior to commencing the approach. The cost 
function used for optimization AS the ITAE which 
accomplishes two objectives. First, it forces the apfroach 
tc be acccmplished in minimum time. Secondly, it insures 
ehat the fuel expenditure will be optimized in the 
eliminaticn cf most overshoot and bang-bang control in the 
dead Zone fortion of the control law. The final value of 
the positicn error must be within the specified dead zone 
and the ternzinal speed must match the reference speed 
(SPDO1). The cost function has the following forn: 


ty 
a= f (te JADX|) dt 


+, 
Table III-4 is a comparison of the optimization -fcuns 
With various initial speeds. The values shown for SW& must 
be multiplied by the speed differential (SPDO2-SPLBO1t) to 
obtain the ccrresponding value of ADX. The max/min values 
show the kand cf values which preduce the optimum cost. 
This cance cf values is attributed to the integration step 
size used in the optimization program. Experience with this 
particular optimization program indicates that errcneous 
values of the switching point are found if the step size is 
nct carefully chosen. The step size may be adequate for 
macegraticn, but net for location of the switching foint. 


The fpeints obtained from the optimization runs are 
Fletted in figure IiI-77. These points define the nenlinear 
Switching curve which must be stored in the computer to 
insure optimal operation of the speed ccntrol for all 
approach speeds. From here there are many frocedure oftions 
Ofen. These cptions have as a goal some usable ferer for 
predicting the optimal switching point for any set of 
initial ccnditions. One may choose linear straight line 
segments with an interpolation routine, cr a closed forn 


Switching curve pclynomial. Due to the availability of a 
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cost function based on 20 min problem time 
all others based on 10 min problem time. 


Table III-4 


Optimization Results 
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hykrid ccnfigured XDS 9300 digital computer and AGT-10 
graphics terninal, the latter course was taken. 


A pclyncmial curve fitting algorithm was used to cbtain 
the reguired pclynomial coefficients of best fit. This was 
done for fpelynomials of order 1 thru 5. The coefficients 
and the sum cf the sguares of deviation frem the original 
Feints are tabulated in table III-5. The selection cf the 
order to fe used is highly deprendent on the degree of 
accuracy required. In the RAS problem, the average error 
introduced for a first order fit is 8.0 feet(1.07? séc), 
while the fifth order fit introduces an average error of 
1.35 feet(C.180 sec). Prior acceptance of errors intrcduced 
by an integraticn (and problem) step size cf 0.8 sec allows 
for use cf a second order fit without any degradaticn of 
Simulaticn accuracy [second order average error is 2.848 
feet (0.38 sec) ]. The graphic display of figure I1II-78 
indicates very little difference in the switching curves for 
second tc fifth crder polynomial fits. Fer the sake of 
accuracy, and owing to the computer control methods of this 
thesis, tke fifth order polynomial fit shown separately in 
figure III-79 is used for determination of the switching 


FCint locaticn. 


A trte test of the control law 1s accomplished when it 
is introduced in a computer program fcr a complete RAS 
Simulaticn. Considering the performance of this ccntroller 
in a complex envircnment of full scale RAS simulaticn allows 


maximum verification of the controller design. 


The scenario for this simulation initially fositions 
the ships such that the ship being controlled starts an 
approach 5 ship lengths (2639 feet) behind the reference 
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ship and displaced 0.4 ship lengths (211 feet) to the right. 
The desired final position is alongside and displaced 0.2 
ship lengths (106 feet). The heading contrel system used is 
developed in section A of this chapter. 


The approach phase 1s accomplished with the _ speed 
desired and speed acquired shown in figure III-80 with the 
corresponding position attainment exhibited in figure 
mer 81. These plots show excellent switching and oftimal 


position attainment. 


The next step is to insure that the fosition keeping 
lecp will maintain the desired position with an induced 
perturbation. This is accomplished by turning the reference 
ship away frem the control ship a total of 15 degrees to 
observe the reaction of the speed control locp. The 
reference ship's turn causes the relative mction between the 
ships to ce altered, making the control ship lag the desired 
position. The nonlinear control system is designed to 
correct this situation as soon as the actual position is 
outside the limits of the dead zone. Figure III-82 displays 
the desired speed and acquired speed for the control ship. 
Figure III-83 indicates that the corresfcending fosition 
deviates frem the desired by 0.0154 ship lengths (8.13 feet) 
at the maximum excursion. This is well within the limits of 


acceptakility for such a drastic perturbaticn. 


The introduction of velocity control was acccmplished 
by combining the simplified engine response of chapter II 
and the speed ccntrol law developed here. By setting the 
speed desired (SPDDES) equal to the output of Function 
SPEREC and scaling the speed error (SPLERR) to the 
nondimensional eguations of motion, the velocity lcop is 
initiated. The auxillary equations added to those presented 
in chapter II are: 
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SPEDES = SPDREC ({ADX,SPDO1,SPD02, Sih) 
CDCT2 = INTGRL(U02,SPDERR*¥LUC) 


These equaticns are introduced in computer program #8 to 
produce figures IfT1I-80 thru III-83. 


Further system study indicates that the reference ship 
speed must be known to a fairly high degree of accuracy. 
Without afpriori kncewlege of the reference ship speed, a 
constant bias is introduced. The amount of bias allcwable 
defines tke fermissible uncertainty in the reference ship's 
initial epfeed. This bias can amount to as much as 0.1 Ship 
lengths (€4.48 feet) for a reference speed inaccuracy cf 2.5 
kncets (0.1 normalized speed). However, it is felt that the 
reference speed in any practical Situation will initially be 
knewn to within 0.5 knots (0.02 normalized speed). fMThis 


mcre practical error will introduce a bias of only 16 feet. 


Cther feedback parameters can be used to offset the 
lack of apricri kncwlege of the reference ship speed. Since 
the reference ship is tracked with a high accuracy rance and 
bearing dévice and the controlled ship's speed is seasured, 
a decoupled multivariable scheme is used tc further refine 
the reference ship speed. With high resciution devices 
presently available{['*], 1t 1s estimated that this can be 
dcne practically to within 0.05 knots (0.002 normalized 
speed). This would bring the offset bias to 1.6 feet; well 
Within previously defined errors introduced by integration 
step size. 
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Thrceughcut the development of the heading contrcl and 
speed contrcl, the scenario has followed the conditicn that 
tke final position would be longitudinally alongside. 
Although this is a good assumption for ships of the same 
type, it dees not account for RAS station differences for 
different ship types. To alleviate this disparity, function 
SPDREC was redesigned to allow pre-planned cffset condition 
to exist. Function SPDOFC cf appendix A is a result of this 


redesign. 


Simulaticn runs, with a change of the speed ccntrol 
function only, resulted in some unstable ccnditions existing 
in the headirg centrol loop. The cause of this pheéencmenon 
stems back to the adaptive gain scheme used and the changes 
made te fcrce the control loop te a steady state value prior 
to aturn. Ey using a favorite ploy of experienced ccnning 
officers, this problem is alleviated. The ploy is to take 
the ship alcngside and then either drop back to station or 
surge forward to station. This method is accomplished by 
setting the initial offset (XOFS) to 0.0. The final desired 
offset (XCFSD) is stored and not used until the ship is 
settled out alongside. It is subsequently used as shcwn in 
the follcwing Fortran code: 


IF (ATIME.GT.450.0) XOFS = XOFSD 


This methcd solved the gain transition problem. it did 
not, however, give a completely stable simulation Cun. 
Unstaktle ccnditions still existed at the end of the turn 
phase. This is not suprising, considering the heading 
ccntrcl optimization method used. The set of gains 


previously fcund were for the alongside scenario only. 


Different interactive forces and moments at the cffset 
position cause these gains to be no longer cftinal. 


By relaxing the control loop in the heading velocity 
feedback gain (VFBG), sub-optinal control at all practical 
orfset pcsiticns 1s achieved.j The gain VFBG was changed 
frem 0.064028 to O.1 in the turn phase adaptive gain 
schedule without significant less cf centrol efficiency for 
alcngside oferaticn (2.3 feet maximum excursion vice 2.0 
feet previously obtained). Subroutine SWICHF of appendix A. 
reflects the gain change and offset calculations required. 
.Ccmputer program #9 incorporates the changes required for 
offset simulation. Table III-6 is a cross. reference Listing 
of the plets obtained. From these figures, the effect of 
different Icngitudinal positicns is readily apparent. An 
offset cf 0.1, eguating to 52.8 feet, causes greater lateral 
excursions when astern (XOFSD = -0.1) of the alongside 
position than when ahead (XOFSD = 0.1). ‘The longitudinal 
position naintainment, however, is essentially the same in 
foe Cases. 
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. 
Note: These plots are the same as those obtained from 
computer program #8 and are not repeated here. 
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Table III-6 
Position Offset Testing Cross Reference 
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An alternative to the eae shown here is again a 
completely acaptive gain scheme which would achieve oftimal 
control instead of the sub-optimal control settled for here. 
The alternative may become even more important if the 
nmcnlinear terms of the equations of moticn are considered. 
This would couple the heading and speed control designs to a 
larger extent than encountered in the interactive forces and 


gctwents. 
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5- Have Effects pmevelccity Control 


The final testing procedure involves validation of - the 
speed ccntrcl system in the presence of waves. This 
perturbaticn testing continues that started in section A. of 
this chapter for heading ccntrol. In chapter II the WX 
force was acdeled thru the intermediate force IF32 as: 


* 


IFi2Z = KCleD2+NC2+KC1¢HX 


By intreducing the force in this way, a severe limitaticn is 
placed cn tke magnitude of the force. In the mariner model 
used, the KC1 coefficient (XDELR) is considered negligable 
cr, at best, cnly 0.00005. This translates, in the original 
equations cf motion, to a maximum speed perturbation of only 
O0.C3255 kts. for the wave amplitude chosen. The second 
drawkack cf this method, with even greater ccnseguences, is 
that the perturbation is introduced before the contrcl loop. 
Delay of the wave perturbation is produced making it out of 
phase with the cther wave force {WY) and moment (WN). 


eee ee ee i ee a ween te ee ee — — _ — 


SPDOES 


Regine bie 102 
Block Diagram of Wave Introduction in Speed Loop 


mimGrdeE eco Dring about uniform introduction cf this 


wave force, its effect is inserted just past the integrator 
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of the speed contrcl loop as shown in figure III-102. This 
is coded in the DSL simulation program as: 


CDOT2 = INTGRL(U02,SPDERR*LUC) +KS2*WX 


A value of ~-1.0 for KS2 will give a maximum wave 
perturbation of 0.85275 kts. (a auch more realistic 
perturbation for the high sea state simulated). Figure 
III-103 portrays the speed desired and speed acquired for 
the apprcach phase in the presence of sea state. From this 
it can be seen that the speed acquired is very dependent 
upon the sea state present. The contrcecl law, however, 
presents a very stable reference for the speed locp which 
gives an approach longitudinal position (DX) plot 
indistinguishakle from that of figure III-81. More 
preminent perturbaticn results are evident in the turn phase 
flets cf figctres III-104 and IIi-105. The speed response of 
figure IIiI-104 allows a maximum longitudinal position 
excursion of 9.5 feet (0.018 normalized position). as 


ccBpared with 8.28€ feet (0.0157 normalized) in calm séa. 
These results show that the speed cecntrol system is 


very stable and corrects well for large external 
perturbaticns. 
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IV. CONCLUSIONS AND RECOMMENDATIONS 


A. CONCLUSICNS 


The results of this design study have been most 
gratifyirg. The Fasic concepts initially ferceived for the 
RAS contrcl have been realized. The decoupled ship ccntrol 
in the RAS environment is a viable and plausible idea. This 
thesis ccntains a workable system for implementaticn of 
ccmuputer ccntrolled RAS. The achievement of 2.3 foot 
maximum excursion for lateral distance while both ships are 
in a tturn, and longitudinally offset by 53 feet isa 
phenomenal achievement. Having this kind of accuracy in RAS 
Ccperations, can vastly increase the safety of this 


ccaplicated and dangerous maneuver. 


The appreach phase of RAS can be a very hair raising 
experience. Night replenishment and sea state ccmplicate 
the "“seargan‘s eve" method now employed in the fleet. Having 
a system that autcmatically handles the apprceach regardless 
of the adverse ccnditions can, again, de nothing more than 


increase the safety of the RAS maneuver. 


Schemes for ccmputer control of nonlinear systems and 
the purpcseful introduction of nonlinear ccntrol laws are 
beccming mere practical with the technological advances in 
Micro precessors. The ever increaSing number of U.S. Navy 
ships with computer systems installed, makes digital 
computer shif control realizakle in the present time frame. 
A gcod wicre ccmputer or an existing installed computer 
(Such as cne used for the NTDS system) can be used in this 
veln. Frocurement of the hardware required for this RAS 
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system can be dissipated over time periods contingent cn the 
funding available. The supply ship requires cnly two 
reflectors fcr the range and bearing devices stationed on 
the recéiving ship. All ships can be outfitted with such 
reflectors at a minimal cost, while the bulk of the hardware 
can be introduced to the ships at regularly scheduled yard 
periods. 


In tke initial conception of this thesis, a secticn was 
rlanned fcr cpen ocean maneuvering. After scme research on 
this facet: of ship control, it was determined that werk in 
this area has already been documented[ ** ]f 25 ]. The 
existence of NTDS outputs for station attainment and single 


ship contrcl systems, made design in this area redundant. 


The concept cf integrated centralized ship contrcel has 
been in the tackgrcund for over a decadef 26]. Although 
given a lcw pricrity due to funding considerations, its 
implementaticn seems to be just arcund the corner 27}. 
However, a review of ref. 27 indicates a lack cf RAS 
capability. Whether this is an oversight in the article or 
neglected in the désign criteria is unknown. If it has been 
neglected in the design, a very real preblem has tfeen 
overlooked. The recent incidences of ship collisions while 
conducting RAS [28] emphasizes the need for inclusicn of 
this very cangerouS maneuver in the “Integrated Eridge 
System." Lack of technology can no longer be used as an 
excuse. This thesis and other research reports [2791] have 
advanced the irzplementation feasibility tc a level that 
cannot be ignored. With these projects finalized in 
practical terms, their incorporation into fleet use is the 


next imperative step. 
A major effort in this area must be made. The ever 


increasing ccmplexities of today's naval shifrs and the leads 


being placed cn the officers and men are such that ccmfuter 
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ccentrol must be used; and used now! We cannot afford the 
luxury of time to prove these system u crth, but must make 
ccncerted effcrts to get them implemented before the lives 
of 300+ men are lost. 


Whenever a complicated system such as a ship in the RAS 
situaticrn is enccuntered, many facets have we be 
concurrently analyzed. This fact has caused inclusion of 
many diagrams in this thesis to illustrate the total 
picture. Each run, with a different ccndition, requires 
Many filets tc analyze the differences in the responses and 
the causes cf the differences. The computer programs shown 
de not reflect the actual run times in the JCL shcewn. As 
Many aS twenty flots were output in these programs in the 
tines listed. Analysis of the actual computation times show 
that the algcrithms run considerably under the time required 
for real time operation. The sampled data rate used in the 
Simulaticns was 0.11 seconds. This is well within the 
realizable data rates available in even the slcwest of 
today's ccmputers and microprocessors. The thrust of this 
consideraticn is that there are nce preblems envisioned in 


ccnverting RAS simulation to real world RAS control. 


BP. RECCMMENTATICNS 


In the heading control design section of this thesis, 
the need for a completely adaptive gain scheme was cited. 
Again in the velccity control section, when a longitudinal 
offset was introduced, this need became even more evident. 
The first and most important recommendation for further 
study is the development of just such an adaptive gain 


scheme. 


The linear equations of motion should rte replaced with 
nonlinear equations to validate the control designs advanced 
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in this thesis. Along with Case the hydredynamic 
coefficients for the Navy's modern ships are required to be 


able to design these control systems for todays vessels. 


It is further recommended that a ccncerted effert be 
made to cktain data on the interactive forces and mcments 
between shirs of dissimilar tyres and sizes. These forces 
and moments gust also be available for sea state conditions. 
In fact, tke whcle area on séa state effects cn the various 
ship types in the RAS situation and in open ccean 
Maneuvering needs attention. Not enough data was available 
fer this researcher to be able to pinpoint sea state effects 
on ship hulls. Since replenishment at sea is rarely 
ccenducted in the sterile condition of calm sea, these 
consideraticns are of utmost importance to allow testing of 


any contrcl system in the Simulation stage of development. 


The intent here is not to imply that the ccntrol 
systems pcertrayed in this thesis are the best for the RAS 
scenaric, but that the procedures used can Lte applied to any 
control scheme desired and bench marked to the ones 
ccrtained here. AS previously mentioned, much meaningful 
research and éesign must be accomplished tc allow system 
reliability and, mcre important, system acceptability Ly the 
Peracers and men who will ultimately trust their lives to 
it. This is a task that must not be taken lightly. 


PAVE 


APPENDIX A 


Due to the lengthy nature of the computer prcgrams 
presented in this thesis, many functions and subroutines 
were develcpred tc simplify their presentation. This 
appendix lists these functions and * subroutines in 
alphaketical crder. The computer programs reference this 
appendix and indicate the placement of the required 


functions and subreutines. 


A Erief description for each listing is given to aid 
the reader in determining their purpose and use. The 
following is a listing of the functions and subroutines 


contained in this appendix in the order presented: 


SCERCCTINE BOXPLX 
FUNCTION DEGRAD 

FUNCTION DELAY 

FUNCTICN FE - RUN A (FEA) 
FUNCTICN FE - RUN E (FEB) 
FUNCTICN FE - RUN C (FEC) 
SCERCUTINE HDGRAS 
FUNCTION KE 

MAIN PROGRAM FOR FUNCTION MINIMIZATIONS (MINIEXPX) 
SUEROUTINE RBMEAS 
FUNCTION RKLDEQ 
SUEROUTINE SLOPES 
FUNCTION SPINIT 

FUNCTION SPDCTR 

PUNCIICN SBDOFC 

FUNCTION SPDREC 

FUNCTION SWCL 

SUERCUTINE SWITCH 


SUERCUTINE SWTCHF 
SGERCUTINE TRANS 
FUNCTION XLIMIT 
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SUBROUTINE EOXPLX 


This sukrcutine was used fer all optimization runs in 
heading ccntrel and speed control. It was programmed 
lecally ard is part of the IBM 360 SSP library at the Naval 
Postgraduate Schocl. A full explaination and description is 
shown in the first few pages of the subroutine listing. 
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FUNCTION DEGRAD 


This function was programmed to ccnvert degrees to 
radians and radians to degrees. A third purpose is to 
ccnvert degrees to a range of 0-360. It is used extensively 
thrcughout tke programs, functions, and subroutines’ listed 


ap this thesis. 
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FUNCTION DELAY 


This function is used as the time delay in the speed 
ccntrcl cptimization runs. It was designed to be used as an 
equivalence to DELY in DSL Simulation.The following 
variables are defined: 


E (1) is the storage array (should be initialized 
before the first functicn call) 


K is the delay step count 
SPDLES is the variable to te delayed for K sters 
P is the flag for delay or no delay 


P>0.0 delay SPDDES 
P<O.0 function cutput equal to SPDDES 


The functicn stores the input value (SPDDES) in E(M) and 
decrements tke value in array E(I) at each call of the 
function until the value is in the positicn of £(1). The 


value is then cutput from the function delayed K intervals. 
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FUNCTIQN DELAY 
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PUNGCTION FE = RUN A 
FEA 


This function is the simulation for heading ccntrcl 
Optimizaticn cf the approach phase. Tt-is. cCalled™ by 
subroutine EFOXFLX. The integration step size is 0.04 with a 
final tine cf 20.0. In this function all initial conditions 
are set to zero except initial geographic locaticn and 
speed. The reference ship maintains a straight course and 
the contrcl ship starts its approach 5 shif lengths astern 
and 0.4 ship lengths laterally displaced to starboard cf the 


reference shir. 


The function is referred to as functicn FEA i1n the 
text. | 


FPR. 


RUN A 
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FUNCTION FE - RUN B 
FEB 


This function is the simulation for heading ccntrol 
Optimizaticn of the turn phase. It is called by subrcutine 
BOXPLX. The integration step size is 0.04 with a final time 
of 20.0. In this function, the following initial conditions 
are non-Zerc: 


contrcl ship rudder angle D2D & ¥(14) = 8.7 degrees 

lateral displacement Y(11) = 0.2 

reference Ship's speed 001 & Y(4) = 1.0 

contrcl ship's speed 002 & CDCT2 = 1.5 (after first 
step beccmes 1.0) 


The reference ship's rudder is. activated tc 5.0 degrees 
between time 4.0 and 5.0. The runs were for port side 
replenishsent. 


The function is referred to as function FEB in the 
text. 


RUN 6B 


FUNC TION FE 
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FUNCTICN FE - RUN C 
HEC 


This function is the simplified simuliation fcr speed 
control cptimizaticn of the switching functicn SW. It is 
called Ey subroutine BOXPLX. The function shown is for 
apEroach speed of 1.1 and a replenishment speed of 1.0. The 
runs were made for various realistic combinations to cbtain 


an cptimur stitching curve. 
The rcur used a step size of 0.04 and a final time of 
10.0. The twce ships were run linearily with only tne 


lengitudinal direction and motion of any ccncern. 


This furction is referred to as functicn FEC in the 
text. 
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SUBROUTINE HDGRAS 


This sukrcutine was programmed to calculate the desired 
heading (FSIDES) for RAS heading control. It uses the 
outputs of sukroutine RBMEAS to calculate this heading with 
gains RSENS and WISENS. The large number of outputs in the 
subroutine call statement were made for ease of DSL printed 
output fcr tracking of simulation accuracy. 


The sukroutine also incorporates a loop to avoid 
computer precision problems in the ARSIN furction. 
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FONCTION KE 


This function is required by all optimization runs. [t 
is the function that contains constraints for subrcutine 
PBOXPLX. No ccnstraints are present, conseguently function 


KE is set egual to 0. 
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MAIN PROGRAM FOR FUNCTION MINIMIZATIONS 
MINIBXPX 


This is a generalized program which calls suwkrcutine 
EFOXELX. It¢ main purpose is input and output of the values 
reguired in the optimization runs. This is referred to as 
MINIEXPX in the text. 
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MAIN PROGRAM FOR FUNCTION MINIMIZATIONS 
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SUBROUTINE RBMEAS 


This sukreutine measures the range and bearing of the 
ferward and after stations which is reguired of sukrcutine 
BCOGRAS. This is done with trigncmetric functions as shown 
in chapter II. The subroutine is specifically designed to 
circumvent any ambiguities usually associated with these 


furcticns. 


It is the basis of the decoupling of the two RAS ships 
that this thesis is based. 
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SUBROUTINE RBMEAS 
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FUNCTION RKLDEQ 


This function is the Runge~Kutta-Gill forth-order 
integraticn used in ail optimization runs. It is prcegrammed 
locally and is part of the IBM 360 SSP likrary. A full 
explaination and description is shown in the first few pages 


of the functicn lJisting. 
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SUBROUTINE SLOPES 


This subroutine contains the table look-up and 
interpolaticrn scheme for the interactive forces and mcments 
presented in the RAS environment. It is long and must be 
pre~compiled fer mest of the DSL Simulation programs shown 


in this thesis. 
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FUNCTION SPINIT 


This furction was designed to aid in initialization 
problems asscciated with the DSL function DELY. The effect 
is that the functicn initializes the delay lcop until it can 
ke self-suppertive. 
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FUNCTION SPINIT 
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FUNCTION SPDCTR 


This function is the speed control used during heading 
control develcpment. It is used directly as the speed of 
the contrel ship with information cf the speed of the two 
ships and the lcngitudinal position ADX. It contains a 
linear furcticn at +1.0 ship lengths to a dead zone of 
40.001 centered about the alongside position (0.0). 
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FUNCTION SPDCTR 
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FUNCTION SPDOFC 


This function is identical to SPDREC except that the 
ability tc offset the alcngside position (0.C) is 
inccrporated. This is the speed ccntrol function in its 
final develcfaent fcorn. 
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FUNCTION SPDOFC 
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FUNCTION SPDREC 


This function 4s similar to SPDCTR except that a 
switching function is incorporated. This is the furction 
used for cptizization of the switching function and 1s used 
in the velccity loop simulated in the velocity ccntrol 
section cf ckapter III. 
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FUNCTION SPDREC 
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FUNCTION SWCL 


This functicn contains the fifth crder felynomial curve 
fit fcr the cptimal switching positicn of the speed ccntrol 
COE. Its range of values for SPDDIF are 0.1 to 1.0 
nermalized speed difference between the two ships. 
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FUNCTION SWCL 
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SUBROUTINE SWITCH 


This stkroutine contains the gains and meéchanisms 
reguired for the adaptive gain schedule developed in this 
thesis. It includes the optimal gains cbhtained frem the 
heading ccntrcl Optimization runs. 
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SUBROUTINE SWTCH 
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SUBROUTINE SWTCHF 


This sukroutine is identical to SWICH except that the 
turn phase gain VFBG is relaxed to allow for cffset 
longitudinal position placement. This is the adaptive gain 
schedule in its final form. 
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SUBROUTINE SWTCHF 
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SUBROUTINE TRANS 


This stbroutine takes the lateral and longitudinal 
geographic displacements and converts’ then to actual 
displacements referenced to the contrel shifts head. This 
is done tc gain a more realistic reference for sukrcutine 
KREMEAS, sukrcutine SLOPES, function SPDCTR, function SEDCFC, 
and function SEDREC. 
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SUBROUTINE TRANS 
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FUNCTION XLIMIT 


This function was developed to allow the LIMIT function 
of DSL tc be incorforated in the optimization runs. [It is a 
Saturation amplifier with a gain of 1.0, and upper limit of 
UL, and a lower limit of XLL. 
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FUNCTION XLIMIT 
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APPENDIX B 


The final fcrm of the simulation program, with all its 
Subroutines and functions, is a very complex and complicated 
maze. fc aid in following its progression, this aprfendix 
ccrnrtains a detailed block diagram of the program with eéach 
variable listed in its computer variable name. Each page 
contains a functional part of the simulation with inputs and 
outputs shecwn cross referenced to their origin and 


destinaticn. 


The fcellowing is a list cf the blceck diagrams ccntained 
in this appendix in the order in which they appear: 


Ship A (Beference Ship) Simulation 
Ship A Heading Simulation 
Ship A Speed Simulation 
Ship B (Ccntrol Ship) Simulation 
Ship B Heading Simulation 
Ship B Speed Simulaticn 
Subrcutine RBMEAS 
Range Measurement 
Eé€éaring Measurement 
Subrcutine HDGRAS 
Heading Cortrol Loop 
Auxiliary Functions 
Yaw Conversion 
Coordinate Conversion 
Feedback Loop 
Rudder Ycdeling 
Ship A Rudder 
Ship B Rudder 
Wave Generator 
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Wave Direction 
Rardcm Variable 
Wave Enccunter 
Wave Components 
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APPENDIX C 


Thrcughcut this thesis subroutine SLOPES has been used 
to output the interactive forces and moments between ships 
in the RAS situaticn. This subroutine, adapfted from ref. 
11, does net ccntain a complete picture of the circumstances 
envisioned. In particular, ship's speeds cther than the 15 
kt. operating point and different ship lengths are not 


accounted for. 


As stated in-chapter II, the speed modificaticn factor 
can easily Le applied for both ships at the same speed and 
otker than 15 kts. with the following expression: 


2 
SPDP = CDOT 


Ships replerishing with different lengths can also be 
inccrporated as shcwn in ref. 1. 


Subrcutine FAMIC listed in this appendix inccrecrates 
these twe ideas along with a better method cf determining 
the interactive forces and moments. The curves of figures 
II-11 and II-1Z were quantized every 50 feet of DX for ail 
the DY curves shown. These pcints were then used in the 
NEGS XDS-¢30C digital computer and AGT-10 graphics terminal 
tc obtain a family of best fit curves. The kest fit 
criteria is tased cn the sum cf the error squared at each 
guantized feint {mcdified somewhat by this researcher's 
e€valuatior cf test fit between pcints to eliminate spikes 
and cther ancmzalies). The results of this curve fit process 
is summarized in tables cC-1 and C-2, which includes 
takulaticn cf the best fit criteria. These pclynomial 


coefficients are based on the DX distance and are ccdéed in 
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Interactive Curve Fit Error Analysis 
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Subroutine FAMIC as YY(1) thru YY(11) and YN(1) thru YN(11). 
An interpclaticn algorithm is used to determine the forces 
and moments at DY points between the curves of each family. 
Although all the computations are based on measurements from 
the ccntrel ship (ship #2), the interactive forces and 


moments are also computed for the reference ship (ship #1). 


Figures C-1 and C-2 are the interactive forces and 
moments cutput to show comparison to figures fII-11 and 
II-12. Tke speed of this run was the operating point of 15 
kts. The ships are of equal length (527.8 feet). 


Linear interpolation of the interactive curves for 
greater than 150 feet DY distance is accomplished from this 
150 foot curve te a vaiue of 0.0 at 200 feet, It therefore 
assumes ne fcrce and moment are present outside the 200 foot 
range. All forces and moments for DY distance of less than 
50 feet are taken as that of the 50 foct curve. These two 
endgpcints are by ne means exact, but will suffice until more 
detailed cata can be gathered. Another inexact endfoint is 
prceduced at the curve families limits of +550 feet. At 
these points, the forces and mcments are forced to 0.0 since 
detailed data outside of these limits was not available. A 
Side effect cf this abrupt truncation will manifest itself 
in the instantanecus commencement of the fcrces and mcements 
during the argproach phase run. The endpoint variations in 
scme of the curves of figures C-1 and C-2 are due in part to 
the curve fitting routine used, but nostly to the 
differences in computer precision. (curve fits were 
calculated cr a 11 digit precision XDS 9300 while the curves 
were flotted cn single precision 7 digit IBM 360/67) 


AS previcusly mentioned, the speed mcdificaticr for 
Other than the operating point cf 15 kts. is cnly completely 
valid fcr tke situation where both ships are at the same 


speeds. Sirce this thesis considered an approach phase 
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where the control ship enters the interactive field at af 
speed guite different than the reference ship, some 
modificaticn cf .- the interactive effects sheula be 
considered. However, exact relationships are not available 


to compute tke required modificaticn factors. 


To dispel any problems with the design of the heading 
ccntrol system, the worst case speed modification factcr was 
chosen. This factor, in effect, considers that the 
interactive fcrces and moments are derived from the ccntrol 
Sid p. This is accomplished in subrcutine FAMIC with the 


fcllcwing fcrtran €xpression: 


SPLE2 = CDOT2Z**2 


As stated in chapter II, it is feit that it is more 
accurate tc consider the interactive forces and moments to 
ke modified Ly the speed of the reference ship, and can be 
coded in svukrcutine FAMIC as: 


SFIF2Z = CDOT1**2 


With the scenaric followed throughout this thesis, this 
expressicr wculd equate to unity throughcut the RAS 
Situation, since the reference ship is maintained at 1.0 


ncrmalized speed (15 kts.). 


For the sake cf error analysis, Simulation of the werst 
case modificéticn is performed. This gives rise tc forces 
and moments 2.25 times what they were in the rest of this 
thesis during a portion of the approach fhase when the 
ncrmalized speed of the control Ship is 1.5. Figures C-3 
and C-4 show the interactive forces and mcments for the 
approach rhase cf the simulation. The comparison plots 
which appear in chapter III as figures III-24 and I1I1-25 
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illustrate the extent of the changes. Most notable is the 
smccther cutput of subroutine FAMIC. This mcre 
realistically portrays the interactive effects in the RAS 
envircnment. Figure cC-5 portrays the geécgraphical flot 
which ccmfpares with figure ITI-26 without speed 
modification. Although differences exist, figure C-5 
illustrates that the interactive effects speed modification 
factor fcr tke worst case does not drastically alter the 
approach phase cutcome. The heading contrcel system cesign 
is still valid in the face of these changes. 


For reference, figures C-6 and C-7 show the interactive 
forces and gcments in the turn phase as calculated by 
sukroutine fFAMIC. Figure cC-8 is the turn phase lateral 
distance frlct preduced. It can be seen from this 
illustration that the maximum excursion is 0.0056 normalized 
distance (2.S6 feet), well within acceptable limits. 


| In sunpmary, the designed control system will accomodate 
even the worst case modificaticn of the interactive effects. 
This insercsitivity to a large range of perturbations, makes 
this contrel system a more viable design fer actual ship 
installaticn. 


3171 


SOOIOG K Poza SAAND oSeyd uANI 


9-D o1NnbTy 
© o0'6 05 °6 no-P 09°C 


(,01 ¥ 99S) BWI 


we *b- 


S2-F- 


~ 
. 


2b 


ST “b- 


CT -t- 


Cv 9 


~ “UOkuum (x T9) 


SZ 


sqZUSUOW N PpeaatTa oAAND sseyud uANL 


L-D ‘eanbt ga 


00 ‘8 03 ‘¢ 


=D 


6 00° | 05° 


(01x 93S) 3wiL 


w 


velo 


O28 


&2°T 


OST 


23 T 


8 °T 


[og 


Z Fowween (x'o) 


els) 


SJOSTIT SOATIOLASRZUL poTJTPOW WorT Aq soueRSTd Tere,zeT eseyd UANI 


8-D) erznbtg7 


OSaU ct 6 65°C CO°G 09 °¢. OC's 


(O91 x 93S) BYIL 


(aq) ADNVASIGA GaHiS32C 


ae 


25 


Ca 
a> 


"2 


$3 


ae 


O-wr<2zuw O> (x'd) 


314 


SUBROUTINE FAMIC 


DQDAOGSCOGOACVODOC YOO COO QGSCOSOSOSCECVBVCOOACOCESGOSGCOG Ceackoet 
mE CO SPUN OP OO OND ION OO FUN GOOD AN COAT INAOPY OA BHAAUOSN OM DADA 
ed ed ed et ted tet IA CIC OIC QUICQUI QI AI CU 09 (9 11 19. 0 HCHO 0) C0) SPSS ES 


ay ey ee ee ee lh eee ee el 
SPE Ea 2 > a a > i > a a a a a 
Sn is We ns Os es os Os a ns Os ss ss sn 
LL UL DR th LL IL LL iUL LLL ih tbh ib i bLibuU th Ub Lh LLih Lh Lt LLU MLR ALLL Lhe) IL LL LL UL LL Uk UL Li 


224DX,ADY,COOT1L, CDOT2ZsY¥1sV¥Y2sYN1sYN2) 


I~ 


© 


CK A RIOR OO IO OI IEG II IO aI a oi aia doak okie ak ds ak ak ofc oka ak ae ok 


DEFINITION OF TERMS: 


| ac 
ths © 
La | ° 
tt Le | 
od om < 
i a) LL. 
alu 
CON fe- 
(a=) fmm] -~t 
Y) J Lu 
<< 
Lt) UL} 
=C > 
— i rm a 
Qu = 
Le 8 | 
CTU 
zz. CG 
Sa? et - 
a 
re a 
ac” Lu 
— mJ 
Tet cy! one) 
YW Wt sash 
UV) <I 
LLoOa a = 
— — () = 1 wm 
Hit. TEE aD] 
WU Zn) a 
LiL © 
me LL UL Li ag 
Z2zZe-OCt: C 
band pag male: 
wZOOAIAS te. 
atl CO WU tt <I 
SE tO ae LLL oom 
LUO OO pret 
AAV LNW mee TT 
—— =O COMM w 
Tor IsTOaMNw VW) 
MMT OL WIL Fi 
ZWINNZooO OW 
LL. UL, #0 /) pene pee CT bo § 
Bao aad = 
mp PC GES e Pare 
CT ke TS TVOOWWO 
aad eek © 40 400 4 a 0 >> Siew) 
OOOwAaIoOaCoOK 
22Zz2RrFfZ22Z2uUus Sse 
Wii O< © 
pe es ee ee | 
QO 
1 IN 
= 
AX > CORINA II 
a wdAQDO OS Kaz 
MM LOO PO 


S15 


SO CRO GIO RR IO ROI IO a IG Io go lok a go x kai OI gO ga ia dak ok iit ae oc ako ake sik ae akc ak ak ak ok 


ic 
eal ee 
=p = b 
© 
<= > 
~ | Oo 
am T< 
<I x< 
OQ ww 
~ «ff => 
cat rome 
te at 
SS -_ 
ae WS) 
a) fol (j 
c= Aa 
wy FT 
(f) tt 
az 
Ss. & a 
— wm A 
MM c+ wd 
Se 2 W 
eave GG. 6 
ao. Wh —t —_ pet C/) 
mam rk WY Ce 
Ez WN Mx 
WwW gf ft 
«ce cow 
ww 9 .-—— 
QO mem AU 
“OZ 
5S a” { ii 
OO * eh OO) 
<{ Le eg on} 
CTC a DMwui 
Cw es COYWU 
Qa LDL je OAiw<t 
A LL. wo 
<{u Onw 
>zZ > mW 
Cel ci © 
ee aw 
a m= O70 LU 
ICD WI W<F QO. 
QwTre~mo Ww 
wm CUI a > 
YW <> OOQ 
WiMrecy) It Uy 


° OLOWOZOO 


U nth Ou <y =” 


NCVE: 


9 
VOVVOO OO’ OO OOO OO OC OQOOQOV OO OO OW 


% 
3¢ 
+ 


CO OO IGOR ROC CK IGG ICO z lok ok IOC IG kk tatoo dol kogodok dtok ied dk kak 4 ak tek 


Ne) 

Cc 

_— 

OO 

Cc 

> | 

8 

© 

wy 

uy 

© 

AN Or 
% eo) 
te HA me Vd 
NN 3 <u 
fe fa OK On 
oo aX 
CQOoa6o *x«KoO 
Ou s« JOd 
aw tt ff +S) 
> 4 uo JC 
ami! a oS 


Hagen J =~ 
Qnaa jt M tl 
see | q LL. 
KIM MY IX ee 


= 


PCLYNOMTAL POWER PRE-COMPUT ATION 


GOO OO GSO OCSCSAdSVOVEOESOBSVWOTSTOMOSCAVDOOOXOE SOeec2ssooecoeockoQ 
PUD SO P= BO OD HO 69 PUNO BD OED OIM FUND PR ODS SQM FUIN OP DOD et GIOTTO ON 
SPs SPP FV LA LAW WI LOA USS LAN NO WO OO WO OO WO OO Be Be fe RPE ME EH DOODDDDDODDOOO? 


Pmt Pret Dat Peet Dat eam Pee Samet ft Peat arnt emad 4 fred ead Pant mat Paes Pt Pee) emt Patent Gad fumed fend Paes ie” badd ae) Pat fell aes eds at ad PM 
SPPSPMWPBSAPTPTPSASASAAAAATSPALTASPLAASTSSATADPASAAPAASZTAS SASL AAT ITST Ss 
Idd {di didd (ete tdi diet did di didi tet i dd efi didi dict dd dit td did died dig etc 
LRU LR LL LR LL LL ER LL ER ALLL UL LL LLL LL LLU EL ULEL ELL EE EL UEL IL UL LL EL ULL LLL UL LLL LLL LLL 


MMO O PLINOM © 
ont ed ed ed ed ed ed dd 
ANI CANT LA SO MH CO OD HH FE 4E 3 Ht He St Ht HK 
HHH HHHRKHHK HH HHH 
3h Se 3 3b He HE IE GE OK OK OK OK OK OK OK OK OK 
> MK OK OK HK OK OK OK OK 
| | 
1 | | | 
OretNMOSSFUVOr © 
FS ONDE SLED FD OS cre pred pred eed pred pad pnd pd ed 
IK HM DK MK DK OK OK DK OK WK OK OK OK DK HK OE OK O*K 


EQUATIGNS OF Y FORCES 


UV 


Nekeoekes| 
aDal 
Qos 
+ © @ 
MOO 
> ae ee | 
te fh 4 
C1>< 
5 3% OK 
+ OI 
woe 
LA + © 
CH LU) + 
Oalw 
Fal 
orca 
LEN OO m4 
ee 0 
Ole 
1 eD 
NO e 


fed wp 
om @ 
Oro 
- @e 
MOS 
< + + 
3¢ [ed 
I< ed 
Ot x 
+ ot iC 
LLiICo =x 
-+- © 
e4tLi | 
fh F WW 
Lh th 
OO 2 
NAILA 
eM) O 
OAL 
1 el 
N® e 


ADEN LTD ts 
e@O e emni em ™ 
Ove D® OWS DVWAO 
+ @el -+ ee - 6 @+ 
MOOW NSO MOON 
Meet Kee Ke Pm 
3b Pert tk Beet it Premio 
wt ets Ket A rt it 
Ox KD OFM Ot MK 
+O + ttt +A © 
Woes + WOO WOr- | 
at OWw D+HO VN+oaw 
COUP ww tlt ee LL ES 
OMuw Ty DetiV Peta 
OMnOD WMA WOmMo 
KOO DOM ODOO 
Mmrore WoO etini@mdco 
2OheH ema etTtna 
Ontos Qty SMO e 
1 eMOQW Il emt + oO 
NO eleHOIO e@ NO e+ 


<x 1 Orrx 1 © 
He O | mts w | 
OX Ox MX O 
QHK atk OF = 
+O] <a + OID 
WiOws QuUuig+ 
LEN fe wt | CT) te peg 
OwOWDOWE 
BOM +#& MO Ort + 
O09 LW et CONT LU 
rtOoroonn 
rt TNO wt? OO 
oN ert 
Oat =Ht OAM CO 
1 eMel ew 
AM AOA - 
x1 of <M] e 
te LAO 46 WES 
AIX +H | 
Ot DOK O34 O 
+ NUK He OU 
tii HF tL: it 
O+-eANtOO+&E 
DwO |} Ole 
SOD + tliest.O + 
xf sf LUM UVLO LU 
OMm A= 0 2 
MOO MB O 
CeO OST enn 
OPTnHOnfH 
- elit + af 
COD eOOQD 
O+- sO0O+ eo 
NTO + OTTO 
ex eA) 


ex +e) 


KX 1Orxx +OTRM + ON 
Se | ttt He Oe tt UP 
AIX © IN KO KAN KE OK 
Ot ett OO + AH Ost KH 
+#$OxXO 1 &# ONIX Ae Ko 
WOx OQWLoOt QOWOt+O 
QtENFtNM+ O11 D444 
MtpOwWlOwOwMu- Ow 
MOAl+ LNDODO + OM +4 
PO WISHES OWI MM Qu 
POM et Oop HO AMO pei 
WO DU tf OF OW Tr 
eOMO o eomMr~ eDOOW 
O Mt HDNOD Ort HO MAN 


1 eN e+ | et el al @ 
MOVNOM AOD OHOnO 
<M | ot etx< | of KX] e! 


HUNOM KH NOM NEM 
QU | wt te QU | IAI | ot 
Ost DKA HOKE HO DK 
wt eS te wt Et OI HH 
moO st 1 WOwe NOt OD 
S + HOw PF FAqQnl +N} 
OW +NnOwO | OWS + 
TO $$ tb PwoO Pe tueN +ui 
tO LL UV SP LU VLA LL OY 
NID tO AM OAM MHtOQM 
NOOPO MD Sines Oal 
eM ofp eDAOm 
ONO eBl-OOOmMmo sr 
+ AN tOFr EOD + SOLA 
MOS 0+A7O4 eOr- ec 
Nl eOvost | oO St ! 


3+0.142822E+02*X1-0.533075F+02%*X2-0.103692E+02*X 340, 3 


2eOOn) 


NPOetesTTo!lSrrol at 


ex+-alx ex 1 


ex 1 AX e 


PH OLA Deal He DHD % Ott 1+ Oa + OHH © 
CON BC QU OK OK OAD CS OKLA) KS OK OOLA AI OK OK SPAS OK ret SP 


Or OFF OXKX ODOHHKHO OF *F OHXHHO 
wFOn FOO tt $A + Nets pete |) + Or | 
wo Woo wood woew woo wool 


ame tn te mr ee SL SF | em OY YH Tm 
et COWL OAS OS UW LL Ot LL ST LU LH © LN LW LL WO ML LU oo 
me TDS TD OS ee FT) ST ed OD WINE IL weer OI OI WO we 
Pe TY PD TE DD HMO Oe ON AIO Of 
Po mt DS = OT et SI) LN 9) HL ODN FOIA D> 


KMN FNM aACQUl 


316 


rtp tN) 


Mt OOS 


DOOVDOCQGOOGOOQAVIACOG©CQQOCSCSS VFO OGOEoScoedTSoo0 200002000 
ICO PUN OP CO NGS ACU PDO PH OD DP OM KN OO PVG PDK Q BANOO SFUVO PH DDS HAM FUO POO 
TKNOADADHOKHVOOOCOVC0O000 Aas et a ata SSAA UNNI NO OAM OC COCAM COO 

met ed ed ed ed eed St edt Sd et et et et ed ted ed eee eS St ed i SS et et ee St St St et tt St St St 
fees Sd met fees Ged Peed eed eng fm rem and famed eee Sed freary acend Pamd ad Sed med Get fed See heel Sang ed Pad a) Dad fe! ed ld PD te Se 
eee ee ek ee eee ee a a a a ae a oe ee ce a ee ne ae oe an in en Pe ae, en an in On, as a, ns a a, es DB, i, Ga, Pn Ds a, 
qeieiei ei dei diief efit didi di cic didi dietitian dedi odie ce 


ag a ee ee a Oe OO LO OO On ON OO OO OO Oe OO Oo Dw 
wnt 


MaOod oMnne eO.f ec a ae eOm e OO ei 0 err om ef C 
eMNlor eMOe OF TO ONN OAM ABM QOM Gr Sean OAIN 
ODm OneAN90 + ce or + 6 + @e| + ee | @ @ f e* {ft ee {ee 
+ee weef MOON MOO MNAOW MOO Eo MOooOo MOAeAO MOO 
MOOoDd nNnqQgown -«—~+ es K+ + >, oe a ee <x +S + <x + > <I (if <1 
Kee Ke te KR K BH he_ thesxX< te Hh RRA RM eR 
thet MAX NAKA NKR Neti NX St Xt aXe KR SX 
Nae AX OF KS OFX Ot KN OF KM OFF KM OF8F K OF KM OFM 
Ort KX OF KQ +4*#0O +17 +OHO tte +t Soe < +tse + eit 
+ mise -HIe+QO Wool Wom woo ! moO WOeO Wee WOa woe 
WOO WOA+ DA+towW M+O0O N+ Ou T+ O40 FTH+O MAtO W+Q 
Oro T+tow Ti+ a Aw i OWF OM Due! NWse Oke MU | MyLts | 
Nude AIF OO NDOWO MOUW DWMoOuo Fou OOWwW BOW MmMmui OGW 
Osttu M~MOWO OOVAe ~OO OOMOag ODOT AO NOOO DHA WAM 
DOmeP ASO UNDO WOW NAOO WAT st OO NTO Mero WOM 
=U DAS SONS Ss HN Of = CO MOON Ms F4UV\CO Ut'cOmD esto OS 
SOO8 aero emo econ ecorhO eof eDO emiN On eN@ 
enw eORF ONO e@ OMO Oore ONS ODN OAM OMA OMM 
Orr BNO e | eo t es [| eDS I eN it eM 1 eo ( et | es 
{ est | e4#QOUONnNO «+ NO 4 NO et NOe NGQe AO se NOo«e NO e 
Oe NOe NOW elaexIior <ISs Kx+O0OT KIOwx<K +0O7K+O0OTR +0 K+ 0 
Io <x~1I1O KIOTxHtO+raq #O+ %Odpet 20! At OF et O+HetO+ %*O+ 
wD I HW I HWIANENKOM NXE NXOK WNKOKNKOXNXKOXNKO NKO 
xO NXO NXKXOxKHaOt AH Dita Or aH OSA EOCaSHOHL SHOOK A Ot 
H$s— Ot+ ae OF SHO +NXO + mt x< - I X< = +N XS +S SX et + NX eet se ot + = >< 
Ns +€NK 4€NKO1UOetO wet WORE wok OWOtOweor OWO+ wet 
Ot WO wWO#FOWO#FEN+S O+0O O + = | O+e latest | Nee ieee OO O+0 
t+Q Ses M4teNehMmwow Awe OwOw PTwOWwdwoWwrwoUrwe Due 
WO st1O OWwaoawdor~+O NM+ Ores AO +F ONT eM 40ND +e NT + 
N- OM+ OM+00RPwN OOw NOwoO rwmetnuih sme ONWw Aa 
Nu OM MOAWRDBTMAMN AQAA COMMON ODONMMO”~rMADODONAM aoOo 
Ks UMD NMOMONOANOQO mshi TONOAH OPODMNNNSBOMADATHOD MOQ 
ON TNS TORO eOsrse ene eDmwown Ort eOvM eMOAO eNO eNO 
own eto cOm0 eOrmnwnd OVO OANWN OAM NO~TNOCAINNOVUe OOM 
NN ONT GONMMNOD I! eA e 1 eM + eNe@ + eO o+ oN 0+ 0H 0+ 0O + eo 
of f e& (| e™ et#HQOnrO —=0M SONG BAONOaAOQNQVAOMO SOM sO 
ey =o. SsOQnaIOrx< |] e| x<f] oe <1 e+ m_—ieli«MjIetxteliM+o K+ e 
< | KX} eb NOM HNOM HUVIM FEMOMNEMNGMENGENENS WO 
hyena Noms WOMXKNK<S 1 st xX [=< NX pas NX | AOU 1 ANC 1 SINK eH ONIX + 
<1 RX LINX LRHOH DK COD ODM OK OK GEAKOQHDKOHAKOUHHD GH ov 
%OXMOH DMOAH OMA + AIX HH + ret o< St te et OC HE eet oS ENC EEN CH HI et OK 
ax tt +N KH +a 3H QIOHO Wotn IUiO tO WE} 3 AWOH SWOH MUG Ht Lc 3¢ 
OF AWOKKWSOHO 1! AH+00 ©0+00 +H Or0OT +200 +H0N+0 n+O 
+$AON ¢€ FOO + SOWVWO+ AWS | DWO+ NOI AWOEI NNO! AWG MuUIe 
WS | PU I MNWOF+QAIr Osetia -—N+FeW OO +r MW +$tbdat+ WOO+ wore OO+ 
V+ (et -UIPO + UiOOasuUi-r DOWN OWUUI ODRONNWINOOUamMesUt Mow 
OWOO OW ODPM NWMOO ON" DOOM Cae PaOPreAMOMOANM—MMD MN 
HOAPOANIDMNDONOOASINN S| arom OMmOM OMMON ANN SS SOMN SND AMON 
IANO SOMMAOTUAO eQNANGA CT Riel Biel Died ePeaiq er ODM eNett erry OM emt CD 
MUS SOR N COR ™DOROO BAMS ADDM OFNOCOss00eeONWN COO 
OtrronrtalOoeth e+ ete + enn + AO + ODS CNS CNM + OO + OM 
oleate ON + est OMOdD ee OOS e OOn e& FOO eON Dan otOd DOO 
Ox e>5FON eOOnm e+ O+ #eOvMr+ 0A N+ eQO0OdD+ OM! eOnmn! *Oote rt e 
- eOD+ @COO+ 2DQ0O0EO TFC | BOYS I! OFS | AOFO+OSOINSTOIOCOTAQ FO 
TOrTTO+rOrOtnt ex | mK ex e—N Ox 1 NK ox tN oX IN ex< IN eX IN ex | 


MEN OX FAN OK HANK OH DR Nt OR 
% OAINH DAO HH Ott NAIK KANN X< 
NX KOK KMONXKOQ OF FO Qt 


Oxext Otte OKHKOV HN | 
tated eee eNat + WOO: 


+m 
Wego 


WOO WOO WOQOOQWHM + +4—0 + | 
De lam + 1-H O+ + NOP CHOW W 
mt LL COD LL ON SP LL LLY OO) et ONO) CO) rete 
COP) B= PLN DS er OPE OO PID WO rT 
DON D> OT SONG rma OT Mw 
TOO STD OV NID SD OPO WO FL COW NI Sm et ot 


ANIM 


ASANO SNS SNOT ANM 


QGUATICNS OF N MOMENTS 


E 


Ph OO wt HONE Ot OO 3H ON HE WHIM ORD CO 
CO) tS OS ONC) ret Soe No et OK ONIO SOK AINSI QQy < 
(ORE HXHO1LOF eH 1 Ot HH LOK (OHH (OX 
+#NO +810 +790 851-9 += 
MIVA Ow UWI HWW OO WUITMO 1WWoe nue 


+ Net | 


Nt 10 OFF O+H Het NGI 


M+ 


mm DS LIU CO ee OO US La 0 UL LU ae 8 Li Uh ee LL et 
tt CO ONO ONO OO SPO MTL OM) OM O 
T(E) ODT OO DN DY ee QD UN eer OU er 1 
a eet NE PPM ZO Met ee D Tet SUN 
P= CVC) COUN D1) ON = De OD LNA OOD Des OO OO OY 


SBOP St NCKQ) etNMO SING SNM) 


SAE 


etN 


(ai Sl olelelelelelelelel ele slolel ele] olelelwlelolelelolel el welol al wl slele) dlalalelulelelelwelel sls) 
QANNOILPFUDDOM OHO ANAM SUDO ODE AN (QPFUVOPM DW HO AAOVF OA OM DDS SAOO SUV 
Ste ST FUL LA VLA WA ALAIN NO WOO OHO ODDO DOPE REE EEE EE DOODDDODO 
SSeqaeqeaseqseqasqqqq eae Hes eat sas aS Sat et St Sst ts Sa eS tS eat Sst sta as St Sa Hs St a et 
PO Oe fet med Ge Ped Pd Pat Pee Omen feat Demet emt pg med Od Gadd ed emt bt md Ded fm © eg Pd Ft my ed fet Fe) eet re en Ped ee ee 
SS. 2eSttStetetsst PASTS SASSP SA SST ASSP AAASAASATDSASASST ARS SATA TASTE 
SW ee eee eee eee eee eee 


ecm. erm eOw od. ef — 
Omar ©O©09MN OFM Bet Orn 
lee !|eet oe | e | @e 
OOo MEO MoCo Moe moo 
<re- <I I xXx+—-eF X< |) <1 I 
He ht Het RMI KH PHK HK 
OxX~ae et XH XH MX MX eK 
OFM OHX DHX OF OKXK © 
-mist +t +t +O +*+O+* 
wood wood woo wo worn 
O+r0O FO OFO OF +O 
UC’ LL + FU + (NLL + OVitr ULL | 
Oe ONW MWNWw Oo NOL 
ODM ~eWHt ~UD NO ern 
OMnN OOM~F Oren NAN OOD 
Anm DOr ND AN NMY 
em at eHDd eO co e2 60 
OMmN Oon OMW Ot Ono 
| ent 1 est 1 eN | e ib eM 
NOe NDe NO aA NO CO eb 
<x+Or~x~ 190 KI OTX<X+ «KIO 
+O + -t5t 0 | %OlnritHzHO #+O! 


N<KOXKNKOQ A XKDQAWN XK OM xO 
OtraetOotst Ott OF AOH et 
+N Xt bP N KM $e OK +H OX 
WUOrOWOSOH WOK QWwOt WO*+ 
asatgedt (| NF Cereal N+ Alor 
AOwOow dow tuew0weEent.& 
mot de et OL | MULE ONAL 1 Or | 
MstWwIteaty SMU MONWSP OU! 
NDOMArP DOCINMmOGAaTONNd —_———— -— 
NOOO N HOO DHO OM ATH uy pt r——t pee bg Oo 
eDOvd enor eOTrest es en <>< tee ee 
OMAnNOnwn QatninNnqnaana TO i i << OQ 
+- eDO e+ Of + ol o+ 0OO+ 0O Pane — i ie | “VU 
HONOQO OW AONOAONAION oO ——— tem ~ NN 
<1 el Xe © KI ol X+- ox<X+ fo -O <x xx ~— 
H+ NOM#H NOMS NOM 4 WOX NO OS) ~~ ee aay 
NX | AN X< +4 Sex J test X< | I< + O HH *%H > > 
Ott OKOKHOKOH OKOH OOH OD o— o_o Ct) 
tm HH IX HH HMI XH HOX + OX =! —— ome O <L<x 
WO * NWO + HWO + AUIO #WO+ aad oe deed | 
D -$OOT + AON +$ HOM FNM + eI — @ ee Be | CAO 
PW! Twa | Dws I CNWOdUO ap > Zz —e = 3 © 
SASSO MNA WOE OWOMN WN tly r ao ee 1) SS CuUVuy 
ONNSAHAOPAWUOAIAU © Finaoo © eo—COC o—_———> > 9 ——_4 
=f LAE AAMAS T OSFNOOINAMY 10 = mt O fan fe - ad | — 
eDDMD OOM #OON ev ein SS o~ ++ —+ + 
ON AAtOMMNDONOFONNOMM - Lene tome tt 200 
- aol + eM + etmetit oO + o- © Cc «dO —_ eee 
NOD 2eDOON! OF eOONNDWN —<I @ eet) ~Z $r-zw Oat 
MWl eOM+ eOM+ eOD! elN+ 8 ~-»-O oO o=- >> >> ow~ 
OFO 1 0-79 +0LT0+0FT00OTFO Se)-F © md pat te Se ee Ase + 
Nex IN ox FN eX $A! ox | aX 4FN eT | ld st + + tt Ie, 
aD TH OAH DHOt DOHOHKH DOH DOA wH LI i > >. 4 ee ee 
Kate XA et IX OK eX | XE ONC f—— © Cmte ee tit Od 
# lL Oe | OHH I OHH OF OXF H DUO aOR Rw— Sew OE ow OF 
=| ¢-HQ +O +eHOl reall eaeN axa ee 8 SS a ed a a a od” e°@ 
T inbed UWOOnwoo wo eae it tee 8 OD De De ef Oe fe ee Oe ee OO 
Pree OVS | Pt eam | =f | 1 UL tt ot Gt Ped te Ee 
ME uw oulle ¢ weenwe: “x sees oe Noles 
OME AQHNQIDONNAHKHNOOD Teer aNnNare ao |] mee renee Soe 
wong weet SOND) er NY ST er SO er OO LN eer et LO er II eee ai —— oe yo are 


MmZOW ST ZO SMS NL ZT SZ SONG 
at AIO) De ot LN Oo PON ST OO KOO et 
Co et NO St) SIC OIC OSES 


COC) LLL te CO eS eS eZ 
<TD OK a OK rat mt et at D> DD Dm De OD De CD Dt ee 


MN (ty wr Fa ‘Oo 


318 


elelelolelelelelalelelwilal sala alae) 
SVHOINUOYTUVOM ODO QaAtMMOas UO 
WOAARKVARKWHAAARANOOOOOOoO 
SN Sete A et mt et et et IO OI OI OU AICI 
pammg head Qered fenad Gered Geng Qed Qe fay frond Gee! fumed Pod fang Gand fem! fey Quad fwd 
SSSASASSA ST STS AATSASAT TSS S 
PO 0 Os Oe Oe 0 ee a ee eo 
LL CLL A 


LAL 
oo LVN 
in | Oo 
WIL t i 
Oo Li ob 
e @ Oo 
ae eed e eo 
3¢ ae 
AUN H% 40 
ee A. ae te 
oo J etjale 
x< x< pa as | a 
C+ 3 xx OD 
ANN % + 
a: 4. =o 
OQQ fo Wo EE em | 
oe. A. OO e 
NY OO. ee 
tt + WY e 
Nl ot HEHEHE SANSA 


Oo eK Z NO mee Kr 

e@¢ Orem JS Free Se ror> 

OO Wil it ax<—a>>> frit 
SZ 


=? ar 
Cie zrnuiuizse Onnvaaenwnnie 
—~ Wt ~D _ 
SNe NCI I eet RAIA 


a2aeiber ce SOK Swe roe S 
Dm Po XL et Oe KO et OO 


@ oO 


319 


COMPUTER PROGRAM #1 


This pregram incorporates the ship dynamics of two 
identical Mariner hulis. These hulls are superimposed in 
space to egllcw fcr comparison of the effects contributed to 
rudder mcdeling differences. [In this particular run a step 


and ramp rudder were compared in chapter II. 


Ancther kenifit of this pregram is tc set up the two 
identical ships required for the RAS simulations in chapter 
III. Basicaliy, verification cf the models in three degrees 


cf freedeaz is acccmplished for the Mariner hull chosen. 


The glcts produced in this run are shown in figures 
imi2 and JIi-:. 


320 


COMPUTER PROGRAM #1 


O@e 

MWY 

— (of: 

Cw 

<<f 

0. a 

==. 

Om 

WW 

2 nl 

Ld LJ 

2a) 
N [a\(an) 
" iis pe | 
ud a a uy 
= ce 
=~ OW ao =) 
= LU uy © 
& = jo ® 
MY | 
i | 
te MN OC 
Se) >> a | 
= < S&S 

aoa ~ = ° 
() > > nw © 
= <I <{ =O fh II 
Tp) a Ma COON ge 

QQ Oo o) 

=< 1 4 OOOuWw 
bons “Oca oO 
ee [2 Aled Oe ex e 
ae Lj LJ iodo e S&S 
= Ft Fe “wnuiwo ul 
= —t py OOWwnN a 
we Cox ->aeamnm~ © 
~— <I <{ NZZANO »- 
wr = re eNO cd 
ww? ae & 8 © 
<< Li Ly WL “MNO ® 
Wu 2 es oe bond (T) ed ot | © 
~ — jae UOWS Il iT} 
uy ——Ooog AN 
~ aus ib @+O *yj”nO 
ad og Loe *Owwx< 
Oo Wwt!iOr+road eae 
~ Lig LL Owe Wt awaee 
OS MY UO>a™ am e 
Ov ae ao xX h—O 
ad OO WO eS anid ii 
NI aa. bat Pe OOO at 
= «86FNY TNMSTOOoO 

Les Ly INA) 0) > 
non Odcac ZOCHO « 
© aS OO! IO 
J Ce KINO wD ett LT © 


ZOOM ad IS iS ee 
mt LJ (oa <— 
med UCU ee RR 
TWN JUIN NNM OfH 
2 AO FrRErFS 2222 UO 
SSNS US ZZ OOOO COCO 
NSO HF tn Ht OOO te 


O e 
—~ NO <q 
oN BE @ ae @) 
<I te ra 
+ =. = 
NI Cas SO ae 
aN UY 
Iq DO —s 
‘ei tonawt 200 
(an SIMON ene O43 4 


A CF AO LK Oi Usd oPO Sees 
2>D> Spear >p>ro> ZY DD OWUI POs 6 eam 
ed O OD KOK Se Oe Ke ZIONS KCK KR OOKXH UD 20O IY YY 
Ce ee ed Da a ae | | | 
OOOH WwW WH MH Han me nn ue 
wil SNS SS SAIN NAO OO) see SOO > Sl 
Tmt) St AINA atest QAI) IMO HWW OB Dee LLL 
YQOaAoaZZaactaodtonodtodtorvyYYvooactoate"Nee 


ATION OF THE COEFFICIENTS 


U 


+ 


SA 


—T* 
oc 
“ee 
<a aw) 
“A! — 
LL LL. WM) 
mee = OO + 4 
-~_+- ~™ == 


—— A -—e OO 
OO Nw ete COO 
QQ <<—t FeO <ter 


=H AN OOOO nd rt et te bee 
NN <f O ot OOOO 
Om— 1 | e OOwwOO 
amt} ot eA © 2 eQniwrx<> 
Oates USIOMOOCasOOm & & 
Ze bm pew bm << (Owes Yt ot 
k-OOO#X #€¥ etd t+ HOO 
HOODOO Sere COLO Sat Kk De 
CDE <S  N HH OD 09 C9 weet em frm wer er 
% tH te ee OR he JOO We 
AMO ww FY OO Le 


mc 

a? Ga? 
NIN awn 
at NJ — 
UL, NO 
mim OO % + 
++ “~\™NM NIN 
NN mom — bP 


ANON QOONANE EE 


NNW << O o— OOOC 
ood | © 0D Oww OO 
1 t &Ne ©2 0QnN2Zx<> 


ONNANS OOOO Orna & & 
Pl ed ed oe SOS a ee esa at 
-OOOk* # Tow +t HOO 
NNUNQAOQDANNN LLYLONA RK Se 
ZO OO <f LD RAI NOD OOCD ww bee be er wee 
aS Ht He ie He ie eM bE OO 
md ed coed ed et mt NY (er et ZOOL 


Stile 6 


foe) 

_" 

Ne) 

aN | 

() e 

—_- eto 
> ale |) a 

CL am se CJ 
Fld EP > SS a, rr 
Se? eee 
Aj J > 
LL wer LI CO < ) 


NAO Ss | OW 
Zt ONQAOA 
OmmN 1 Were 
ANOS NICO #€ 
We 0D >= Let 
WON *<fDss 


Axe encase 


mt met OOD UE met C9 0909 et ee HE KPO eee OO ON SS IO =—(OAWM I oe 
SOOO Sse Hh tne RO eK VY ODOONNAN HH He eon NN SS Oh On I 
Wop td ee ett tet SS LU ee meee NINN ZN ZZ WH NLU OL oe EN OU ae SI 
matt CO Ce am am et bm ae a et et Se OTN EEE CO 0 et em ee et I OW OO LO ®OOO 
QNHaetetQHOQAOCOVOCOW WVSACLTAHANANNNOOQOOAOCOCONWOoOOoONW NW SFeelrv~ort iwnetess 
LUMO DOQOAOD SFO Ost RBLLLHNMOODOADOOANOQANN St ZD DPDeanwLZzrvaeat 
rat bot parg etd <I OD) (_) <I C) _) CD) DK Sw KD em TL) /) nt et pet ed oe eet <TC) <0 CO ae OKO ENN LO QIN > > 


S22 


DYNAMI 


C 

=) 

S 

o 

a) 

a) 

= 

< 

~ 

at 

COO 

MW e 

OOo 

o if A come 

IY) = Ajom 

Od YO 

Sw acctQ~ 

JQ 2 EeZzO 

~F* NNO 

wp © a> O 

00 aN ee 

et @ LU LULL LU 

So COO 

ei Cee 

te DDE 

OY NNT O 

Ud eee oe J 

IaQO SNS CW 

> & cnn ee 

ow @ ete AIC ~~ 

=m OOOO cw 

me jf Sszese O 

Se 

I= QOQ0Q0 wUW 

> 

FoJjwdadad/ 

et a ed el 

ALlOdeef<ct ic 

© QOOQU ato 
Ww jk = 
are — 
Q. Zhe wad = 
> Pa mm 
Ivox Lu 
Ma Oa. — 


OT-SYSIN DD * 


6.0 


8.0 


0.0 


-~6.0 iO 8.0 6.0 


100.0 


0.0 


328 
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¢.0 


INSERT TWO /* CARDS HERE 


COMPUTER PROGRAM #2 


This pxrcgram models a practical rudder response fcr a 
Mariner shir type. The rudder limits (stops) are set at +30 
degrees anc the rate of response is limited to +2 
degrees/seccnd. A scale factor (LUC) is introduced to 
modify the respcnse to match real time of the mariner hull 


chosen. 


Twelve passes thru the program are accomplished to 
cconfcrm to different sets of initial conditions and final 
desired rudcer conditions. The plots produced in this run 


are shown in figures II-5 and II-6. 
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COMPUTER PROGRAM #3 


This pregram models a reduced order (first order) gas 
turbine froprision plant for an input-output relatianship. 
The program does not scale the plant to the mariner hull 
used. This was done when introduced into the main 


Simulation prcegram first listed as computer frogram #8. 


The time delay (P) is assisted in initialization by a 
dual feed into the system; one thru the delay itself and one 
directly intc SPDIN. The programs can be modified to ccmpare 
a family cf curves. by introducing the following sequence 
into the TERMINAL region: 

INTIGCER NUMB 

IF (NCUR.EQ.NUMB) CALL ENDEW(NPLOT) 

IF (NCUR.NE.NUME) CALL RERUN 

NCUR = NCUR + 1 
where NUME is the number of curves desired (less than or 
egual tc 10) which is set with a PARAM statement. The 
comparison is dore on the conditions set in the terminal 
region [i.e. decrement or increment the system gains (€g. G 
pa + 0.CZ) }j. 


The plet procuced by this rcun is shcwn as part of 
figure II~10. 
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CCMHPUTER FROGRAM #4 


This pregram models a simplified wave Simulation 
ccmpesed cf two superimposed sinuscids (fundamental and 
seccnd harmcnic) and a small random wave. The aodel is 
inherently scaled - to the Mariner nondimensional 
characteristics. Introduction of these waves is accomplished 
in computer progran #7. Subroutine DEGRAD is shcwn in 


appendix A. 
The sea state force plots in the dimensions cf the 


three degrees cf freedcm produced by this run is shown in 
figures [I-17 thru II-22. 
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COMPUTER PROGRAN #5 


This ficgram uses the mariner hull model facse 
introduced in computer program #1 and the control system 
designed in chapter fII to simulate the afprroach phase of 
BS . The subroutines and functions that are to be inserted 
frcm appendix A can also be dore in object ccde by changing 
the word FORTRAN to OBJECT and placing pre-compiled decks in 
the same locations. In fact, due to the lcng length of 
subroutine SIOPES, this must be done to be able to run the 
Sipulation with the DSL default job control language (JCL) 
presently irstalled at the Naval Postgraduate Schocl IBM 
SoC 767. 


The flicts produced by this run are shown in figures 
TiII-7 thru i1I-13. By changing the gains and introducing 
the foilcwirg code, the plots of figures III-14 thru III-19 
are prodtced: 

CDILES = S.0*STEP (8.0)-5.0*STEP (9.0) 
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COMPUTER PROGRAM #5 
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COMPUTER PROGRAM #6 


This pregram combines the approach and turn phases of 
ccuputer efpregran #5. The added subroutine is a result of 


Simulaticn requirements to switch between adaptive gains. 


This rrr froduced the plots of figures III-22 thru 
IlII-34. By substituting the initial conditions of table 
IIi-3, tktis program produced the plcts cf figures III-35 
thru I[IIT-€4. 
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COMPUTER FROGRAM #7 


This pregram combines the calm sea simulaticn of 
ccmputer fregram #€ with the wave simulaticn of computer 
pregram #4 tc Simulate the model and control system in a sea 
state. The waves are introduced thru the rudder 
nendimensicnalized coefficients as shown in chapter II. 


fhe flets produced are shown in figures ITI-66 thru 


III-73. Figure {11-65 was produced with the same program by 
setting wWl=1.5. 
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COMPUTER PROGRAM #8 


This pregram incorporates a fifth crder pcelynomial 
curve fit speed ccntrol switching functicn te give optimal 
lengitudinal fositioning. The scenario is the same that was 
used in tke cesign of the heading control development. The 
low order secdel of the gas turbine propulsion flant was 
used. 


The flects fproduced are shown in figures III-8C thru 
1-83. 
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COMPUTER PROGRAM #9 


This pregram introduces a longitudinal pesiticn cffset 
capability. The sethod takes the control ship tc the 
alcengside esition until 450 seconds into the run. After 
that time, with the ship steadied, the cffset position 
desired (XCFSD) 1s switched to the desired offset. This 
method negates some of the transient oscillations which 
cause urstakle conditions in the apfprcach phase. A 
seccndary change is the use of subroutine SWTECHF instead of 
Sukroutine SkTCH developed in the heading control section. 
This new subrceutine relaxes the heading velocity feedback 
gain (VFEG) tc allew turn stability in the turn phase. 


The plots produced by this pregram are shown in figures 
mer—oe¢ thru Iiz-101. 
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COMPUTER FROGRAM #10 


This pregram incorporates the sea state first 
prcgrammed in computer programs #4 and #7. The WX wave 
force, hcwever, is introduced at the end of the propulsion 
locp: to allew more realistic perturbations. This is the 
final fcrm cf the complete heading and speed ccentrol 


systems. Tc cun this without a sea state, set WFMA tc 0.0. 


This pregram produced the plots shewn in figures 
esos thru [Ii-105. 
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